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SUMMARY 


The method treats wing configurations with trailing vorticity sheets of 
arbitrary shape. Induced angles normal to the wake are analytically formulated. 
These, factored by the vorticity in the sheet, form the induced drag integral. 

The vorticity is represented as a Fourier series in loading, then the unknown 
Fourier coefficients can be determined after specifying the constraint conditions 
then minimizing with respect to the loading Fourier coefficients. Exact analy- 
tical solutions in terms of induced drag influence coefficients can be attained 
which define the spanwise loading with minimized induced drag, subject to speci- 
fied constraint conditions, for any nonplanar wing shape or number of wings. 

Closed exact solutions have been obtained for the planar wing with the constraints 
of lift plus wing bending moment about a given wing span station. Compared with 
a wing with elliptic loading having the same lift and same wing bending moment 
about a span station, the induced drag can be of the order 15 percent less for 
the wing with the inboard minimized-sol ution loading. Example applications of 
the theory are made to a biplane, a wing in ground effect, a cruciform wing, 
a V-wing, and a planar-wing winglet. For minimal induced drag the spanwise 
loading, relative to elliptic, is outboard for the biplane and is inboard for 
'the wing in ground effect and for the planar-wing winglet. The theory can be 
applied to determine the loading for minimal induced drag of wings in formation 
flying, banked wings, flying in or near wake, and linked wing tips. No-roll 
can be an additional constraint. A spinoff of the biplane solution provides 
mathematically exact equations for downwash and sidewash about a planar vorticity 
sheet having an arbitrary loading distribution. 


INTRODUCTION 


Induced drag minimization theory is in a rapid development and use phase. 
This is because of widespread interest in performance and structural advances 
that possibly may be realized with unconventional aircraft configurations, or 
by imposing structural and/or performance constraint conditions. Numerical 
solutions of minimization theory have been programmed for computers, which 
extends the application range of the theory. Unconventional aircraft include 
such configurations as a wing with winglets; a biplane; a tandem wing; a cruci- 
form; a V-wing, swept box wing; a strutted wing; a skewed wing; and a nonplanar 
wing designed to lessen wing and fuselage bending moments. A structural con- 
straint condition can be that the wing bending moment about a given span station 
be minimum. A performance constraint condition can be an aircraft flying near 
ground, a surface flying vehicle, or an aircraft in a flying formation. 

Minimum induced drag solutions for planar wings with different structural 
constraint conditions have been developed by several investigators. Dr. L. 
Prandtl (ref. 1) in 1933 obtained a solution using the constraint that the 
integral of the local bending moment be specified. Compared with an 18.1 
percent smaller span wing with elliptic loading, and the same lift and constraint 
magnitude, the induced drag is 11.1 percent less, in reference 2, the constraint 
condition is the wing-root bending moment. The solution for spanwise loading is 
obtained after the induced angle function is determined. Compared with a 25 
percent smaller span wing with elliptic loading, and same lift and wing-root 
bending moment, the induced drag can be 15.6 percent less. In reference 3 the 
constraint condition also is the wing-root bending moment and the solution is 
somewhat similar to that of reference 2. For approximating wing structural 
weight, the integrals of wing shear force and of wing local bending moment are 
used as constraint conditions in reference 4. With these two constraints, the 
induced drag can be 7.1 percent less when compared to the 13.8 percent smaller 
span wing with elliptic loading and same lift and constraint magnitude. 

Numerical solution vortex-lattice methods for computing the optimum span- 
wise loading which results in minimum induced drag, have been developed in 
references 5 and 6. These are quite adaptable to complicated nonplanar 
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configurations. In reference 5, the optimum loading calculation is based on 
Munk's theorem III in which conditions are defined on the induced velocity. 

In reference 6, the numerical procedure is generalized in terms of constraint 
conditions. An application of a nonplanar wing numerical method is made in 
reference 7 to the investigation of wing-winglet characteristics. Examples of 
nonplanar wing analytical or classical theory solutions for loading are given 
in references 8 and 9. 

The principal objective of the present study is to develop a generalized 
analytical induced drag minimization theory to which aerodynamic and configura- 
tion constraint conditions can be imposed. The term generalized implies direct 
solutions are possible for arbitrary configurations without having to resort 
to complex conformal transformations. Here generalized also means that mini- 
mization is done without resort to other induced drag theorems such as Munk’s. 
Aerodynamic constraint condition examples are, for a given lift, for a given 
rolling moment, and/or for a given wing bending moment about an arbitrary span 
station. Configuration constraints include, multi-wing, wing-winglet, or 
flying near wake of other aircraft. It was found that a generalized theory could 
be formulated by expressing the spanwise loading distribution as a Fourier 
series in which the Fourier loading coefficients are unknowns. With constraints 
imposed and minimization, equations are evolved in terms of the unknown load- 
ing coefficients and induced drag influence coefficients. Solutions of these 
equations leads to an analytically exact evaluation of the loading coefficients. 
In an age of numerical solutions, available exact solutions provide an accuracy 
check. The present theory also serves as such a check for the evolving numeri- 
cal minimization methods. 


SYMBOLS 


A aspect ratio, b 2 /S 

a n , a n * Fourier coefficient of spanwise loading, equation (1) 

b wing span 

C. ratio of wing bending moment coefficient to lift coefficient, 

C mbb /C L’ ec ^ uation 01) 
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'be 


' D i 


'mbb 


av 


D i 


nr.* 


f l » f n 


nn* 


'rrw 


bb 

m,n,n*,N 


C b of wing with elliptic loading, f-j/47r 
induced drag coefficient, D^/qS, equation (13) 

lift coefficient, L/qS, equation (3) 

rolling moment coefficient, L/qSb, equation (78) 


wing bending moment coefficient about lateral point n h5 M hh /qSb, 
equation (7) 

wing chord 

average or mean wing chord, S/b 
section lift coefficient, equation (5) 
induced drag 

induced drag influence coefficient for planar-wing winglet, 
equation (164) 

induced drag efficiency factor, equation (15) 

induced drag efficiency factor of rolling wing, equation (79) 


wing bending moment influence coefficients for planar wing, 
equations (9) and (10) 

dimensionless circulation, r/b V, equation (1) 

downwash influence coefficient for biplane, equations (67) and (118) 


induced drag influence coefficient for biplane, equation (70) 


sidewash influence coefficient for biplane, equations (119) and 
(B12) 

planar wing bending moment constraint parameter of C, and e, 
equations (26) and (37) 


planar wing bending moment constraint spanwise loading function, 
equation (31) 


planar wing bending moment constraint parameter of r\ , equation 
(33) CP 


wing lift, also, rolling moment 

induced normal velocity influence coefficient for wing winglet, 
equations (165) and (166) 

wing bending moment about span station 

integers in Fourier series for loading and in influence coefficients, 
odd only for symmetric loading 
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induced velocity influence coefficient for cruciform wing, 
equation (125) 

induced drag influence coefficient for cruciform wing, equation 
(128) 

induced velocity influence coefficient for V-wing, equation (141) 
induced drag influence coefficient for V-wing, equation (143) 
dynamic pressure,^- pV 2 
wing area 

spanwise coordinate along nonplanar wing surface 
spanwise semi span along nonplanar wing surface 

wing bending moment constants, equation (171) 

arbitrary constant, proportional to how much wing bending moment 
constraint is imposed 

limit condition for t, equation (29) 

free stream velocity 

induced sidewash, tangent to surface 

induced downwash, normal to surface 
lateral coordinate, from midspan station 
vertical coordinate 
induced angle normal to the wake 

circulation, also dihedral angle 

angle between wings in cruciform, also, dihedral cant angle of 
wi ngl et 

influence coefficient constants used for evaluating I n and I ,, 
equations (B39) and (B40) !u 1 

dimensionless vertical coordinate, 2Z/b or Z/s e 

dimensionless lateral coordinate, 2y/b or y/s g 

span station about which bending moment is taken, equation (6) 

span station at which winglet starts, dq = cost^g 

semi spanwise center of pressure, equals 4C^ when ^ = -rr/2 


spanwise trigonometric coordinate, cos -1 n 


Subscripts 

a v 

b 

c 

m,n,n* 

s 


average 

bending moment 

wing with elliptic loading 

series integers 

symbolizes, from vorticity sheet 
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MINIMIZED INDUCED DRAG OF PLANAR WINGS WITH THE 
CONSTRAINTS OF LIFT PLUS WING BENDING MOMENT 


Relations between spanwise loading quantities can be established in the 
wake or Trefftz plane of the aircraft. Then solutions reduce to two-dimensional 
problems and are independent of the chordwise shape of the wing and of how the 
loading is developed. Quantities dependent on spanwise loading include, induced 
velocity normal to the wake, induced drag, spanwise center of pressure, wing 
bending moment, and rolling moment when the loading is antisymmetric. When 
constraint conditions are imposed and induced drag is minimized, then from the 
relation between loading and induced drag, the optimum spanwise loading can be 
evaluated. A simple method for establishing the relationships between the 
spanwise loading quantities is by expressing the spanwise loading as a Fourier 
series in which the Fourier loading coefficients are unknowns (see eq. 1). This 
method permits a manageable, viable, analytical solution, regardless of the 
complexity of the configuration or of the constraint conditions. The method 
can be demonstrated by showing the solution for a planar wing with simple 
constraints. With the constraint that lift is specified, then (by eq. 15) 

3 ^d . 

C °i = ^ (a ! 2 + n = Z 3n2) ’ ^7 = ^ am>m>1 = ° s G = a l Sin * 

that is, the loading for minimum induced drag for a given lift is elliptic 
loading. Also, for a planar wing with the constraint that rolling moment is 
specified, then 

3C 

C D • = a 2 2 + f " a n 2 )’ aa^ = m a m )m^2 ° s G = a 2 sin 2(J) 

i n= I f 2 m 

which is the loading for minimum induced drag for a given rolling moment. 

In this section solutions are developed for the constraint conditions of 
lift and wing bending moment about a given span station such as at a wing strut 
connection. It can be noted as pointed out in reference 2 that the induced drag 
of a wing having a given lift and a given spanwise load distribution is not 
affected by the compressibility of the air at subsonic speeds. At supersonic 
speeds an additional drag associated with the formation of waves arises. 
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Force and Moment Coefficients in Terms of 
Fourier Coefficients 


The wing spanwise loading distribution expressed as a Fourier series has 
proven to be a versatile analytical tool (see e.g. refs. 10, 11, 12, and 13). 
In Fourier series 

c c 

G(<|>) = bv = 2B~ = I a n sin n# (1) 

where for symmetric loading n is an odd integer, and <j> is related to span 
station by 

n = 2y/b = cos <f> (2) 

Lift coefficient is 


c - M 

L qS qS J 


b/2 


pVrdy = A jGdn = A 

-b/2 -1 


r 

G sin 
o 


4> d(j) 


(3) 


then with equation (1), equation (3) becomes 

r - TrA n 
2 ® i 


(4) 


When C^ is known the loading distribution can be written 


~ G = - (sin <f> + t sin n<j> ) 

W| b| 7T ^ ^ d J 


(5) 


'Lav 


n = 3 
odd 


Wing bending moment coefficient about the span station 

% = 2 Vb = cos *b 


(6) 


IS 


M 


bb _ 1 


'mbb qSb qSb 


r b/2 


(y - y b )pvrdy = j 




(n - n b )0dn 


(cos (p - cos <p 


G sin <j> d4> 


(7) 
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then inserting equation (1) and performing the integration 


'mbb 8 


(fi + f n f -f) 
n = 3 n ai 
odd 


where 


fi = 2(sin <j> b - ^ sin 3 ( 


’j-j cos <f> b ) 


( 8 ) 


(9) 


f = 1 r sin (n-2)<J>h sin (n+2)4>h 2 cos 4>h sin (n-1 )4>h . 

n n L n-2 " n+2 ~ n-1 

2 cos <J>b sin (n+D^hi. sin (n-2)^h 2sin n , sin (n+2)4>h /-, n \ 

n+1 J" Tn-^yfrTTJn ' (n-l)n(n+l) n(n+1)(n+2j (10J 


At 



n-1 

-4(-l) 2 

n(n 2 -4 ) 


(10a) 


An important parameter is the ratio of wing bending moment to wing lift, 
the ratio of equations (8) and (4) 


c b = 


^ . 1 ( f + f „ f JL ) 

C L n = 3 n ai 

odd 


From 


01 ) 


where f 1 and f are given in equations (9) and (10). When the loading extends 

over the semi span, then <f> b = tt/ 2, and is 1/4 of the spanwise center of 

pressure, n rn - Then, combining equations (10a) and (11) 
cp 


rl 


n-1 


'cp 


c ft c j _ 4 fi o % 

C c r) ° r| 3 tt l 1 - 3 2 n2-4 

u L c av n = 3 n 4 

odd 


, -.s 2 a 
(-1 ) _n -i 


Induced drag coefficient is 


c 

D. qS qS 


rb/2 

pVr dy 

A 

2 

f 1 

G “w dn = f 

-b/2 

• 

-1 J 


G a , sin 
w 


$ d 4> 


( 12 ) 


(13) 
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where a w is the downwash angle in the wake given by 


f 77 


a - — 
W TT 


dG 


J o 


d<j>* d(f>* 1 «? . 

bos^cos» = iW n * f 1 n a n* Sln "*♦ 


Inserting equations (1) and (14) into equation (13) leads to 


(14) 


r a C, 2 a 2 C, 2 

C D. =54. E n a 2 = — j— E n(J-) = -V 

l 4 n = i n "A n = 1 'a x ' irAe 


(15) 


where for symmetrical span loading 
1 


e = 


~ co ,3n, 2 

1 + e n (~) 
n = 3 31 

odd 


(16) 


Solution for a /ai Fourier Spanwise Loading Coefficients 

The problem is to evaluate the infinite number of a n / a i Fourier coefficients 
such that the spanwise loading is defined which has minimum induced drag per wing 
bending moment about n^. Minimal combined induced drag and wing bending moment 
per unit lift can be realized by minimizing the parameter C^/e. From maxima 
and minima theory for a function of many variables, the necessary conditions for 
an extreme value are that all the partial derivatives with respect to each 
variable be zero. Taking partial derivatives of C b /e with respect to a 3 /a ls 
a 5 /a ls .... a m /a ls results in 


3(C b /e) 

3(a m /ai) 



3 C, 


3e 


3(a m /ai) e 2 afa^/aj) 


= 0 


The partial derivatives of equations (11) and (16) are 


(17) 
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( 18 ) 


ac. 

b _ m r: 

a ( a m/a^ ~ 47r m 


9e 

3 ( a m / a i ) 


2m( a m / a i) 


= -2 e z » - 1 

a i 


3^2 2 

[1 + " n (jS-) ] 


n - 3 
odd 


With equation (18), equation (17) becomes 

a 


m 


m 


7i — — f + 2 m C, = 0 
4ire m b a x 


which can be solved for a m /a 1 . Then the solution for the Fourier loading coef- 
ficient, with m referred to the general n, becomes 


n 

a i 


8-rreC, 


09) 


Equation (19) is in the form of a constant times (-f ). Let t be an arbitrary 
number, then the minimization condition is 


— = - t f 

n 


( 20 ) 


A better parameter for the minimization process is C^ 2 /e since it indicates 
induced-drag bending-moment efficiency. However, the minimization of the term 
C^ r /e where r is any constant leads to (a n /ai) - -(r/87TeC b )f n , that is, a 
constant times f and can be expressed as in equation (20). 


Wing Bending Moment Minimized Solution Loading Characteristics 

Solution in infinite series . - Equation (20) inserted into equations (5), 
(11), (12), (14), and (16) leads to the following spanwise loading characteris- 
tics : 
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( 21 ) 


C o C ?A 4 , w 

p— — = p — G = — (sin 4> - t z f sin n<j>) 
u L c av l L 71 n = 3 n 

odd 


cp 



+ 3t 


n-1 


f (-D 2 

n = 3 n2 - 4 
odd 



( 22 ) 




t 

sincji 


z 

n = 3 
odd 


n f 

n 


sin n<j>) 


(23) 


r _ mbb 1 , , . 

C b - = (f i ‘ kt) 


ttAC^ 1 + kt 2 
i 

where 

k = z n f 2 

o n 

n - 3 

odd 

where f is given in equation (10) and in equation (9). 


(24) 

(25) 


(26) 


A limit on the magnitude of t. - A limit can be defined by the condition 
that the spanwise loading remains positive at all span stations. This condition 
is satisfied by requiring that the slope of the loading distribution be zero 
at the wingtip. Taking the derivative of equation (21) with respect to <j>, settin 
to zero, and solving for t at <f> = 0 


t 1 “ 


Z n f 


n = 3 
odd 


(27) 
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where t is the value of t at which the loading is positive at all span stations. 
If t exceeds ti the spanwise loading becomes negative in the region near the 
wing tip. The summation in equation (27) is readily found by expanding n f . 

From equation (1(D) 


? n f = sin 4> b + J sin 3 4*5 + g- sin 5 <f>^ + sin 7 4> b + 

1 . _ 1.-7 

- ■g sin 5 ^ - j sin 7 4^ - 
- 2 cos 4> b (^ sin 2 4> b + | sin 4 <j> b + ^ sin 6 4> b + 


n = 3 
odd 


- ^ sin 4 4> b - -g- sin 6 


. ) 


n 

= sin 4> b + j sin 3 ^ - sin 2 4> b cos <|> b = y sin 


(28) 


then equation (27) becomes 

t = 3 = 3 

1 2sin 3 <j* b Z ^_ nb 2)3/z 


(29) 


Solution in closed functions. - The infinite series summations of equations 
(22) and (26) converge quite rapidly, however, those of equations (21) and (23) 
do not, especially that of equation (23). Summation methods are developed and 
presented in appendix A which includes the functional evaluation of the sum- 
mations given in equations (21), (22), (23), and (26). Using equations (A60), 
(A74), (A87), and (A95) the closed function form of equations (21) through (26) 
become 


c * c 


C L c av 




n 2 ) k (1 - k Q t) 


(30) 
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where 




“»' r ‘ & - ! » - v > V! 


(31) 


cp 3n 


tr 0 


- kit) 


(32) 


where 


k i - ?r (1 " f % 2 ) (1 - n b 2 )^ + \ n b 4 cosh" 1 ~ 

b 


(33) 


ttA 

2C l a w 


- 1 + t 


H 


o, | n I ^ n b 


, 7T ( I n I _r i b ) » I n I >n b 


(34) 


r — mbb _ 1 / _p i i \ 
c b - Tf - *T (fl - kt) 


(35) 


where f : = j (2 + n b 2 ) (1 - n b 2 )^ - 2n b cos -1 n, 


k - (1 - 6n b 2 + 3n b 4 + 2n b 6 - 6n b 4 In n b 2 ) 


(36) 

(37) 


e = 


1 


C 2 
L L 

ttACj-j 1 +kt 2 


(38) 


The functional behavior as n b approaches unity of the coefficients k, ki, 
and fi can be made as done in the development of equation (A75). 
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For n b 1 : 


k = ^ (1 - n b ) 4 [1 - fo (3 + n b ) (1 - % )] 
ki = fg (1 - n b ) 5 / 2 0 + % ) H [1 - f 0 - n b )] 

fi ' if * <! - %> 5/2 tl * Iff 0 - - b >] 


Solutions in closed functions for n b = 0- - In this case the solution 
simpl ifies. The condition that bending moment is taken about the wing root is 
that = 0- Then the k, k l9 and f x factors simplify considerably. For = 0; 
k = 2/9, ki = 1/6, fi = 4/3. With these values the aerodynamic characteristics 
are, for = 0 


= - [(1 + I-) O - n 2 ) h - t n z cosh'i 1 ] (40) 

C L c av * 3 |n| 

-cp-fcn-l) HD 

3^% = 1 + 1 ( f - w|ri! ] (42) 

c b = h n - 1> (43 > 


Numerical values of loading characteristics for various nb and t. - In this 
wing bending moment minimized solution several generalized constants appear. 

These include k, k l9 and f x . The term k Q is a generalized function of the span- 
wise coordinate. Values of k, k ls f 1# and the ratio k/f x have been computed 
from equations (37), (33), and (36) for a range of values. These are presented 
in table 1. Values of k Q determined from equation (31) appear in table 2. 
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The spanwise loading distribution is given in equation (21) in the form of 
an infinite series, and in equation (30) in the form of a closed function. It 
is of interest to know the magnitude of the loading harmonics, that is, the 
factors of sin 3 4 >, sin 5cj>, . . . Values of f computed from equations (10) and 
(10a) are presented in table 3. These data show that the optimized loading is 
made up primarily of elliptic loading and the sin 3<f> harmonic, and much less 
due to higher harmonics. 


The value for t such that the spanwise loading gradient is zero at the wing 
tip is ti given by equation (29). For these t = t : values, the aerodynamic 
characteristics are listed in table 4. These are computed by use of equations 
(32), (34), (35), (38), and table 1. The parameter (C^ 2 /e)/ (C^/e^ , c denotes 
t = 0, is a measure of the induced drag efficiency compared to a wing with 
elliptic loading having the bending moment ratio C ^ c /C ^ . The optimized loading 
is elliptic when t = 0. The parameter of the induced angle in the wake has a 
constant value from n = 0 to then it is a straight line joining the points 
at n b and at 1 . 

The bending moment ratio is obtained from equation (35). Then 



(45) 


Solving for t 

1 = W7 (1 - c^ ) (46) 

For comparison of data it is convenient to keep the parameter C^/C^ constant. 
Then t is evaluated from equation (46). The aerodynamic characteristics for 
C bc /C b = 1-2, and 4/3 are presented in tables 5, 6, and 7. These values 

are determined with the t of equation (46) and using equations (32), (34), (35), 
(38) , and table 1 . 
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The spanwise loading distribution due to the specified conditions on 
loading in tables 4 through 7 is obtained by using equation (30) with the listed 
k Q values in table 2, and t values of tables 4 through 7. Values of spanwise 
loading are presented in table 8. For positive values of t the loading is 
inboard of elliptic loading. 

For -> 1 , mathematical limits need to be taken. With equations (39) and 

(29) the terms kjt, fjt, kt, and kt 2 can be evaluated and are all zero at n b = 1 • 

Then n = 4/3ir, and e = 1 , that is the values for elliptic loading. However, 

C P _ r, 

(kt/fi) -> 5/8, and tt(1 - n^t *> 3 tt t/ 2/8 ( 1 - n b ) 2 > as n b -* 1 • Thus for the t of 

equation (29) as n b -* 1 


ttA 

2C l aw 


(n) 1 , 


ttA 

2C l a W 


(!)->!- 


Since from equation (39), (k/fi) 
(46) , as n b + 1 , becomes 


3 tt /2 C bc 8 

8(1^7* ’ Cb "" 3 

+ ( 5^/6)(l - n b ) 3/ \ 


(47) 

then t of equation 


t 




3/2 

5(l-n b ) 3 /2 



(48) 


where C b /C bc is a specified value. With the t of equation (48) then, as before, 
rj C p 4/3tt, e 1 , and « w (l) at n b = 1 approaches minus infinity. The values 
for n b = 1 are incorporated into tables 1 through 7. 


Comparison with Elliptic Loading 

A measure of the efficiency of a wing with inboard spanwise loading can be 
attained by comparing this wing with a wing having elliptic loading and the 
same lift and wing bending moment. Since elliptic loading is the most efficient 
in terms of induced drag it serves as a good standard for comparison with the 
inboard loading characteristics. The induced drag ratio and ratio of aspect 
ratios are 




2 





(49) 
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where subscript c refers to the wing with elliptic loading. The wing bending 
moment coefficient ratio is 


'be 


. W C L 

^7 


M 


bb 


L„ b 
c c 


M bbc L 


(50) 


for the same lift and bending moment, and since for elliptic loading e c = 1 , 
then equations (49) and (50) reduce to 



(51) 


D. 

i 


D ic 


= 1 ( jL _)-2 
e ^b ' 
c 


(52) 


The ratio of aspect ratios depends on S/S . It can be noted that when these 
wing surfaces are the same then friction drag will remain about equal. In 
equations (51) and (52), C b and e are given in equations (35) and (38), 
respectively. The wing bending moment parameter for elliptic loading is 
equation (35) with t = 0. Then 


C 


be 


fl{n bc ) 


4tt 


(53) 


where n^ c = ^bc^^/^) 1S the span station of the wing with elliptic loading at 
which wing bending moments are taken. Two conditions on can be made. These 
two moment conditions are pictured as follows: 
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wingtip, for wingtip, for 
elliptic load- inboard loading 

ir ^ \ 
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(54) 


Then for the same y b 
b 

n bc b c n b 

For these two cases of moment definition, equation (53) becomes 

fi(n b ) 

Same C bc = — ^ 


(55) 


r 

Same y b : C bc = 


(56) 


with equations (55), (56), and (35), equation (51) becomes 

C . b = 1 

Same n b - 57 k_ 

c 1 - f 


(57) 


Same y b : r - f _ v kt “ 


(58) 


Solving for t leads to 


Same n b : t - ^ (1 - F }^) 


(59) 


Same 


v 


t - 


1 




k/fi k b/b 


(60) 


where k and fi are given in equations (36) and (37), and in table 1, and 
f i (bn b / b c ) is 

b \_ 2 ro ,/b r- / b \9-1h n b — -wb 

T n h' ' T L 2 + 

c 


fl ({j- n b ) = | C 2 + (f- n b ) 2 ] [1 - n b )‘]* - 2 n b cos-^ n b ) (61) 

c c c c c 
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When b/b c is specified then t can be determined from either of equations 
(59) or (60) for a given value of n b - With this value of t, induced drag 
efficiency, e, is computed from equation (38), induced drag ratio from equation 
(52), n C p from equation (32), C b from equation (35), and induced velocity is the 
wake from 


ttA 

2C l a w 



w 


wc 


= 1 + t 



hl<n b 

(1^1 -n b ) 9 |n|>n b . 


(62) 


When t is specified, such as for the condition of wingtip zero-slope 
loading of equation (29), then b/b c can be determined from either of equations 
(57) or (58). An iteration for b/b c is required in equation (58) when t is 
specified. 

The data presented in tables 4 through 8 applies also to the case of same 
n b when comparing this data with that due to elliptic loading. This can be 
seen by noting that equations (46) and (59) are identical when C bc /C b = b/b c- 
From the relations of equations (51), (52), and (62), the data of tables 4 
through 8 can be written as 




ttA 

; 2C l a w 
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wc 


C b Ve 


= 1 (— — ) ~2 = 

e K b> 
c 


D ic 


(63) 


which relations apply when comparing with elliptic loading which has the same 
value of v n bc = V 

Data for the case of same n b is tabulated as a function of b/b c in table 
9. A similar table for the case of same y b is presented in table 10 with 
corresponding spanwise loading values given in table 11. 
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MINIMIZED INDUCED DRAG OF NONPLANAR WINGS 
WITH EITHER THE CONSTRAINT OF LIFT OR OF LIFT 
PLUS BENDING MOMENT 


In the yz-plane or cross section plane the trailing vorticity system is 
assumed to analytically follow the downstream projection of the aft view geo- 
metry of the aircraft. Thus a staggered wing biplane will appear as two 
planar vorticity sheets separated by a vertical distance. The term nonplanar 
means that the vertical position of the vorticity sheet or wake is a function 
of the lateral coordinate. Nonplanar can also mean multiple independent 
vorticity systems which, however, can influence or interfere with each other. 

A flat straight vorticity sheet is the simple planar wing. 

The general approach to a minimization solution is first to analytically 
determine the induced angles normal to the wake. These angles factored by 
the vorticity in the sheet form the induced drag integral. The vorticity in 
the wake is represented as a Fourier series in loading. Then the drag integral 
develops into an infinite series in terms of unknown Fourier loading coefficients 
and induced drag influence coefficients. Constraint conditions are specified 
on the Fourier loading coefficients. Then the induced drag series is minimized 
with respect to the Fourier loading coefficients, that is, partial derivatives 
are taken with respect to a given a . Solutions for the Fourier loading coef- 
ficients, a n , are obtained from these minimization equations. Constraint 
conditions can be as simple as that for a given lift, then a-| is constant in 
the derivative process, or for a given rolling moment, then a^ is constant, or 
for wing bending moment with a given lift, then the induced drag series times 
the bending moment series in a n is minimized. Thus, exact analytical solutions 
in terms of induced drag influence coefficients can be attained which provide 
the spanwise loading for minimized induced drag, subject to specified constraint 
conditions, for any nonplanar wing shape or number of wings. 

Biplane and Wing in Ground-Effect Solutions 

The biplane configuration provides an example for a nonplanar minimization 
solution. The wing in ground-effect is the same but with opposite sign of the 
vorticity from the lower wing. Shown on the right are the trailing vorticity 
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sheets from wings a and b. The normal induced 
downwash at sheet a, due to both trailing vor- 
ticity sheets is 


w(y,Z) 


sheet a 


— - 
sheet b 0 



t”dr a (y r z) 


dr b (y ] ,Z) 


w(y,z) 



1 



dr. 

d V: 


dyi 


y-yi 


1_ 

2n 


b b /2 
• * b b /2 


dr b 

(y-Yl) 3yY dyi 
(y-yi) 2 + Z2 


(64) 


For the same wing span and vorticity strength, in trignometric coordinates, 
equation (64) becomes 


a U) = d. 
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7T 

G'(^)d^ 

COS^i-COScJ) 


1 
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r 

7T 

(cos^t - coscj)) G 1 ( i ) d4> i 

(cOScf^-COSc})) 2 + C 2 
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(65) 


where G and cos<f> =n are defined in equations (1) and (2), and 5 = 2Z/b. The 
Fourier series for G from equation (1) is 


G(<J>i ) 


co . - , \ 00 

L a sin n<f> L , G ' ( 4> x ) = 2 n a cos n^ 

n = 1 n n = 1 n 


The first integral is 


1 

IT 


7T 

cos njj djj _ sin n<f> 

cos<j) i-cos<() sin <p 

Jo 


and the second integral is defined by 


I = - 
n t r 


^(cosf 1 -njcos n^i d<J>i 
(cosfi-n ) 2 +' 

0 


( 66 ) 


(67) 


where n~ cos $ and I is evaluated in equations (B5) through (Bll) in appendix 
B. With equations ( 66 ) and (67) equation (65) can be reduced to 


a w 


n = 1 


( s4n_ji^ + j ) a 

v sin <f> n' n 


( 68 ) 
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The induced drag coefficient (eq. 13) becomes 


= A 
' D i " 2 


A 


G a w sin<j)d<}) - 2 


l a sin n+ S + V " 

n = 1 n n*= 1 Sln + n 


( 69 ) 


A / TT 


• i r\ i Jl w o w w , 

sin<j)d<j) = «■ («■ e n a 2 + E a E n* a 
^ ^ n = 1 n n = 1 n n*= 1 1 


I * sinc^sin n<j>d(f>) 


Define 


nn* it 


1 ^* sincfisin n<|> dcf> 


(70) 


then with equation (4) for C, , equation (69) becomes 


ttA C Di 
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■? 





n* 

1 




(71) 


These lift and drag coefficients and aspect ratio are based on wing area and 
wing span of one wing only in this equal winged biplane. The I nn *' s represent 
induced drag influence coefficients due to a second wake vorticity sheet. The 
I * integral of equation (70) is evaluated in equations (B25), (B26), (B27), 
(B39), and (B40) of appendix B, also in table 12 for 5 = 1/2. 

For symmetric spanwise loading only odd numbered n and n* apply, then 
equation (71) can be written 



Equation (72) applies for a biplane with equal wings. Since I nn * denotes the 
influence of one wing on the other, then for a monoplane wing in ground effect 
the I's in equation (72) are all preceded by a negative sign. For elliptic 


23 



loading, a n = 0 for n>l , then e" 1 = 1 + 21 for the wing in a biplane, and 
e” 1 = 1 - 21 -ji for the wing at h/b = c/4 above the ground. In ground effect 
the ground vortex images are in opposite sense to the vortices of the model. 

a n /a-) coefficients for minimized induced drag. - Objective is to minimize 
equation (72) for a given lift but without any other constraint. Following 
the procedure leading to equation (17), taking partial derivatives of 1/e 
of equation (72) with respect to a^/a-p a s/ a i 5 • • •» a m /a-j > resL ^ts ln 


= 2 m — + 2(m I, + I , ) + 2 Z ml ^ + 2 z n* I 
a m /a l a l lm ml n = 3 nm a l n*= 3 


nn‘ 


odd 


odd 


— = 2[m I, + I , + m — + z (ml + n I ) — ] 
a-j L lm ml a-| n = ^ nm mn' a-| 


(73) 


odd 


where m is an odd integer greater than one, and the I's are from equation (70) 
or appendix B with corresponding subscript integers. For minimal induced drag, 
equation (73) is zero, then 


m It + I t + m 
lm ml a 


m 


!<-,.*» W a7 = 0; m=3 - 5 > 7 - • 

n = 3 1 

odd 


.N (74) 


where N is an odd integer made large enough to insure accuracy. There are 
( N- 1 ) / 2 m-equations with (N-l)/2 a /a-j unknowns. A simultaneous solution of 
these equations leads to evaluation of the a n /a i ratios. As an example, for 
N=5, equation (74) appears as 


3 1 1 3 + I 31 + 30 + 2I 33 ) ^ (3I 53 + 5I 35 ) ^ = 0 

di 3 

51 1 5 + I 5 l + OI 53 + 5 I 35^ if + 50 + 2 I 55 ) 4 = ° 


> 


(75) 


For the wing in ground effect a negative sign must be inserted before the 
I values of equations (74) and (75). Examination of equation (74) shows that 
minimum induced drag is not attained with elliptic loading. For elliptic 
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loading, a n = 0 for n>l , then equation (74) is not satisfied because m I^ m + 

I 1 is not zero, except when r, - 00 where 1 = 0. With the a n /a-| ratios known, 
then T/e can be determined from equation (72), induced angle in the wake from 
equation (68), wing bending moment ratio from equation (11), n rn from equation 
(12), and spanwise loading by equation (5). 

an/a] coefficients for mi nimized induced drag and wing bending moment. - 
The ratio C^/e is minimized as in the development leading to equation (17). 
Then 


3(C b /e) _ r a ( 1 /e ) 


ac, 


e d ( a / a -, ) 
' m 1 


= 0 


(76) 


where a(l/e)/a(a /a-j) is given in equation (73) and 9C t/ 9 ( a n/ a l^ is glven in 
equation (18). With these values, equation (76) becomes 


a N a n 

ml. + I . + m z (m I nm + n I ) — + m t f 
1m ml a-j _ ^ mn a-| m 

odd 


= 0; 


m=l ,3,5, . 


.N 


(77) 


where t is a constant similar to that in equation (20), and f is given in 
equation (10). Equation (77) is equal to equation (74) plus the added term 
m t f . The solution procedure for determining a n /a-| is the same as that 
of equation (74). 


Antisymmetric loading a p /a? coefficients for minimized induced drag 
The wing rolling moment coefficient is given by 
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c * ■ 


qSb 


Gridn - ^ 


-1 
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G sin 2<j>d<f> = g— a^ 


(78) 


where L is the wing rolling moment, and G is defined in equation (1). With 
this C and relation, equation (69) can be written as 
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( 79 ) 
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For antisymmetric loading only even integer values of n and n* apply. Then 


] = 1+21 + — 

e r 1 dL 22 2 


n = 4 
even 
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n = 4 
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^ nI 2n + 2I n2> J 


n = 4 
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E n*I 
n*= 4 
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nn 



(80) 


where e^ is a drag efficiency factor for loading due to rolling. The partial 
derivative of equation (80) leads to the equations for minimized induced drag. 
Thus, 


aO/eJ a N a 

—7 — t — y = ml 0 + 21 0 + m — + z (ml + nl ) — = 0, 
a^/^) 2m m2 a 9 n = A v nm mn' a 0 


m = 4,6,8 . . 


2 n = 4 
even 


(81) 


which applies to the biplane. For wing in-ground effect the signs of the I's, 
are changed. 

Example numerical solution with c= 1/2 of biplane and wing ground effect 
model s . - For N = 9, equation (74) is written as 

a-, a r a- 


3I 13 +I 31 +3 ^ 1+2I 33^^ 3I 53 +5I 35^iy f ^ 3I 73 +7I 37^aj + ^ 3I 93 +9I 39^a^ 
5I 1 5 +I 51 + ^ 3I 53 +5I 35^aT +5 ^ +2I 55^aT + ^ 5I 75 +7I 57^aT + ^ 5I 95 +9I 59^aT 


i -| u-j 

7I 1 7+ I 71 + ( 3I 73+ 7I 37)^+( 5I 75+ 7I 57)^+ 7 ( 1 + 2I 7 7 )^+( 7I 97+ 9I 79)~- 


9I 19 +I 91 + ^ 3I 93 +9I 39^aT f ^ 5I 95 +9I 59^aT + ^ 7I 97 +9I 79^af f9 ^ +2I 99^a! 


= 0 
= 0 
= 0 
= 0 


> (82) 


l l 


1 


where I nn * is given in equation (70) and in appendix B. Using the I nn * for 
C = 1/2 given in table 12, equation (82) becomes 
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a 3 a 5 a 7 a 9 

139103 + 3.348175 — - .155810 — + .009726 — - .004216 — = 0 

a i a i a i a i 

004910 - .155810 — + 5.179748 — - .078668 + .010699 = 0 

a l a l a l a l 

a 3 a 5 a 7 a 9 

001489 + .009726 — - .078668 ~~ + 7.077076 ~ - .044194 zr = 0 

3-j 3- j 3-j 3-j 

003196 - .004216 ^ + .010699 ~ - .044194 ~ + 9.008733 = 0 

3-j 3-j 3-j 3 -| 


> ( 83 ) 


J 


The simultaneous solution of these four linear equations provides the Fourier 
loading coefficients. Convergence accuracy can be assessed by also doing the 
solution of the top three, top two, and top equations. These are 


4 eq. solution 
3 eq. solution 
2 eq. solution 
1 eq. solution 

From equation 



.041649 .002196 -.000241 

.041649 .002197 -.000243 

.041648 .002201 

.041546 

(5), the spanwise loading distribution is 


a i 

.000370 


(84) 


— = — (si n<(. + S sin n<f> ) (85) 

C L c av " n = 3 a l 

odd 

where a n / a i values are given in equation (84), and 4> = cos" 1 !!. This loading 
is that for minimum induced drag of a biplane wing with the wings vertically 
separated by c = 1/2. This optimized loading is outboard relative to elliptic 
loading. 

With the c = 1/2, I nr) * values of table 12, the drag efficiency factor of 
equation (72) becomes 
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] ± .422074 + l 
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± [-.278206 ~~ - .009820 — + .002979 
a. a. 


9-^ 9^ 9^j 9^ 3 7 9 q 9 

-L + -r- (.348175 ~ - .274717 — + .019452 - .008431 — ) + — 

a^ a^ a^ a, a^ 

(.179748 - .157336 — + .021397 ~) + — (.077076 — - .088388 

a.| a.] a^ a-j a-j 


a q 2 a q 

-^) + .008733 (— — ) - .006391 — ] (86) 

a i a i a i 


where in the ± signs, plus applies for the biplane a^/a^ 's and minus for the 
ground effect a^/a-j 1 s. 

The spanwise center of pressure of equation (12) simplifies to 

^ i a r n 3' 


'cp 


_ (1 + 3 !l . 1 !5 + 1_ __ 
3ir u 5 a ] 7 a, 15 a. 
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With the a n /a-| values of equation (84) used in equations (86) and (87), 
the 4 through 1 equation solutions compare as follows: 



4 eq. 

3 eq. 

2 eq. 

1 eq. 

el 1 i ptic 


sol . 

sol . 

sol . 

sol . 

loading 

1 

e 

1 .416271 

1 .416272 

1 .416273 

1 .416294 

1.422074 {88) 

n cp 

.434873 

.434879 

.434885 

.434993 

.424413 


which shows 0.4080 percent less induced drag than that with elliptic loading. 

The minimum induced drag and root bending moment solution for the biplane 
is obtained from the simultaneous solution of the equations of equation (77). 
That is, to add mtf to the left side of equations (82) and (83). For the 
wing root case, = 0, and using equation (10a) 

m=3 . m=5 d m=7 . m=9 d 

mtf m “ 5* " FT 1 45 1 1 = ’ TT* 
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Choosing t = 1, these values are added to equation (83), and its solution 

for a / a n results in 
n I 

a M n=3 n=5 n=7 n=9 

-3- = -.195876 = .031633 = -.012113 = .005932 (90) 

an 


and for odd n > 9 ^n 4(-l) 2 / 

a-j n(n 2 -4) 

where for n > 9, a n /a-| is obtained from equation (20). With equation (91) 
the following summation is: 
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n - 11 a l 
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J71 4 n = 11 a l 
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Then with the a^/a-j values of equation (90) and t = 1 


(93) 


£ = 1.613214, 


— = 1 .134410, 
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n = .372175 
cp 


= 1.422074, n cpc = .424413 


r n K _ 0 
c b b n 


- -5E- = .876764 
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(95) 


'be 'epe 

With equations (49) and (50), for same lift and same wing root bending moment 

*ic 


jp = 1.140558, ~ = .872038 


(96) 
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The spanwise loading distribution for these minimized conditions is equation 
(5) in which the a^a-j's are those given in equations (90) and (91). These 
results show the minimum induced drag for a given root bending moment due to 
t = 1 inboard loading of a wing in a biplane with c = 1/2. Compared with 
the biplane wing with elliptic loading and same lift and wing root bending 
moment, it has 12.8 percent less induced drag, 14 percent greater wing span, 
and spanwise center of pressure is 12.3 percent less, that is loading is inboard. 

For ground effect solutions the I nn *' s ln equation (82) change sign. This 
sign change can be taken into account in equation (83) by replacing the 
numbers 3.348175, 5.179748, 7.077076, and 9.008733 by -2.651825, -4.820252, 
-6.922924, and -8.991268, respecti vely. With these numbers substituted into 
equation (83), the 4 through 1 equation, simultaneous solutions result in 




a 3 

a 5 

a 7 

a 9 



a l 

a l 

a l 

a l 

4 eq. 

solution 

-.052494 

.000675 

.000136 

-.000331 

3 eq. 

solution 

-.052495 

.000676 

.000134 


2 eq. 

solution 

-.052495 

.000678 



1 eq. 

solution 

-.052455 





(97) 


where as before in equation (84), excellent convergence is demonstrated. The 
a n /a- 1 's of equation (97) inserted into equation (85) gives the optimized span- 
wise loading distribution for minimum induced drag of a wing in ground effect 
at h/b = c/4 = 1/8 above ground. This optimized loading is inboard relative 
to elliptic loading. For elliptic loading, a n / a-| = 0 for n > 1. 

With the a n /a-| values of equation (97) used in equations (86) and (87), 
the 4 through 1 equation solutions compare as follows: 



4 eq. 

3 eq. 

2 eq. 

1 eq. 

ell iptic 


sol . 

sol . 

sol . 

sol . 

loading 

1 

e 

.570628 

.570629 

.570629 

.570630 

.577927 


.411014 

.411008 

.411004 

.411056 

.424413 
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which shows 1.2627 percent less induced drag than that with elliptic loading. 

The minimum induced drag and root bending moment solution for the wing 
in ground effect is obtained by the simultaneous solution of the equations of 
equation (83) which is changed by substituting the four numbers listed above 
equation (97), and by adding equation (89), t = 1, values with changed signs. 
The solution results in 

a n=3 n=5 n=7 n=9 

- 11 = -.357148 = .050279 = -.013734 = .005717 (99) 

a i 

n-1 

and for odd n > 9 a n 4(-l) 2 /-,™\ 

17 “ n(n 2 -4) (100) 

These values are inserted into equation (93) for t = 1, then 


1 = .826212, n = .330000 1 

e cp 

[ 001 ) 
±-=. 577927, r 1cpc = .424413 J 

= 1.429613, V° \$_ = _ 777544 (102) 

C bc “ n cpc 

With equations (49) and (50), for same lift and same wing root bending moment 
h D. 

~ = 1.286101 , = .864307 (103) 

i c 

The spanwise loading distribution for these minimized conditions is equation 
(5) in which the a n /a-j 's are those given in equations (99) and (100). These 
results show the minimum induced drag for a given root bending moment due to 
t = 1 inboard loading of a wing in ground effect at h/b = c/4 = 1/8 height. 
Compared with the wing in ground effect with elliptic loading and same lift 
and wing root bending moment, it has 13.6% less induced drag, 28.6% greater 
wing span, and spanwise center of pressure is 22.2% less. 
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The antisymmetric loading solution for minimum induced drag for the 
biplane with c = 1/2 is obtained with equation (81) in which the I's are the 
even n and n* integer values of I nn * listed in table 12. Then for an N = 6 
solution 


-.158416 + 4.256811 — - .107350 — = 0 

a 2 a 2 

.004993 - .107350 — + 6.121520 — = 0 

^2 


(104) 


Then 2-equation and 1-equation solutions are 


2 eq. sol . 
1 eq. sol. 


= .037211, 


— = .037215 


= -.000163 


(105) 


The induced drag efficiency parameter of equation (80) with even integer 

values of n and n* for I * of table 12, becomes 

nn* 


— = 1 ± .215435 + 2(^-) 2 + 3(— ) 2 + [-'.158416 — + .004993 ~ + 
6 a 0 a r\ d a d a 

r 2 2 2 2 


.128406 Cr-) 2 + .060760 fe-f - . 1 07350 ” ~] (106) 

□ 2 q 2 u 2 ^2 

where in the ± term, the + sign is for biplane, and the - sign is for ground 
effect values. Then with equation (105) for 

2 eq. sol., — = 1.212487; 1 eq. sol., -= 1.212487 (107) 

e r e r 

that is, the same to the sixth decimal place. This is 0.253% less than the 
elliptic loading equivalent sin 2 <j> loading value which is 1/e = 1.215435. 

The spanwise loading is given in equation (1) which for antisymmetric loading 
is 



oo n 

z — sin n<f> 

n = 2 a 2 
even 


(108) 
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where from equation (78), a£ = 8 C^/ttA. The spanwise loading resulting from 
inserting a n /a£ values of equation (105) into equation (108), is outboard of 
the sin 2cf> loading. This minimized induced drag loading is for a biplane with 
only antisymmetric loading. 

The antisymmetric loading for minimum induced drag for the wing in ground 
effect at h/b = d 4 = 1/8 height from ground is obtained with equation (81) 
but with sign changes on the I's from table 12. 


.158416 + 3.743189 — + .107350 ^- = 0 

a 6 

-.004993 + .107350 — - + 5.878480 -- = 0 

”2 z 


(109) 


2 eq. sol . 

1 eq. sol. ~ = -.042321 

a 2 

Using equation (106) 


a 4 a 6 

— = -.042368, — = .001623 

a 2 a 2 


(no) 


2 eq. sol., — = .781205; 1 eq. sol.,— = .781213 

e r e r 


( 111 ) 


which is 0.428% less than sin 2$ loading value which is 1/e = .784565. The 

spanwise loading is equation (108) with the a n /a 2 values of equation (110). 
This loading is inboard of the sin 2<f> loading. 


Unsymmetric spanwise loading is a combination of symmetric and antisymmetric 
loading, that is, it includes both odd and even integer values of n and n* in 
the loading series. A solution for minimized induced drag starts with equation 
(69), written as 


C = — ( ? 
D i 4 n - 1 


" V + 2 
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n = 1 


n*= 1 


* ! nn* a n*> 
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Then 
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D, 


air 1 = f [m a m + n ? , (mI nm + nI mn } a n ] ’ m = 3 ’ 4 ’ 5 • • • 
m n = 1 


(113) 


For a minimized induced drag solution, equation (113) is equal to zero, then 
the equations are solved simultaneously to obtain a^, a^, . . . a n< The coef- 
ficients a-j and a^ are specified for given values of C L and of equation 
(4) and (78). For a single wing, the I's of equation (113) are zero, then 
a^ = 0 for n > 2, which means that minimal induced drag is with a-j sin $ and 
a^ sin 2(j) loading combinations. 


Flow-Field Solution of a Flat Vorticity Sheet 

A spinoff of the biplane drag minimization solution provides mathematically 
exact equations for downwash and sidewash about a vorticity sheet having arbi- 
trary loading distribution. The downwash at point (y,z) due to the vorticity 
sheet is the second integral of equation (64), then 


w s (y, z ) = ^ 


b/Z (y-yi) &-dyi 


\y-y i) 2 + z 2 


(114) 


J-b/2 


where subscript s symbolizes, from sheet. Similarly, the sidewash is the 
horizontal component of the induced velocity, that is, z/(y-y^) times the 
integrant of equation (114). Thus 


v s (y> z ) ■ 27 


b/2 



7 -b/2 


(115) 


where the v^ is positive in the y-direction. In terms of the dimensionless 
coordinates of equations (1) and (2), equations (114) and (115) become 
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y— ( n ® c ) 
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n 7 t 
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(117) 


fTT 


(n.c) = -c ? 

n = 1 


n a„ - 

n it 


cos njh d<h 

(C0S<j>!-Tl)2 + \' Z 


where a^ are the Fourier coefficients of loading, given in equation (1). The 
integrals of equations (116) and (117) are defined in equations (67) and (B12) 
of appendix B, then 



n I a 
n n 


v 

s 


V 



n J a 
n n 


(118) 

(119) 


where for arbitrary values of n and c, I n is determined from equations (B5) 
through (BIT), from equations (BIT), ( B 1 3 ) , (B14), and (B18) through (B20). 
It should be noted that a sign change in n influences the sign on I and 
as shown in equations (B6) and (B16). 

For n = 1 the loading distribution is elliptic and equations (118) and 
(14) can be expressed as a ratio so that 
w 

S _ t f i 


where w w is the downwash at the wake and is constant spanwise for elliptic 
loading. The elliptic loading solution for w s / w w is given in the literature, 
for example, solutions by conformal transformation in references 10 and 11, 
however, it appears as a function in elliptical coordinates. The general 
solutions for w /V and V s /V such as those of equations (118) and (119) were 
not found in the literature. Computed values of equation (120) for various 
n and c correlate with values in the table on page 150 of reference 11. 


Cruciform Wing Solution 

The cruciform wing or x-wing is a nonplanar type configuration. An aft 
view of this wing and coordinate system is shown on the next page. The induced 
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velocity normal to the upper right wing 
panel due to the crossing wing with dr 
vorticity is 

dw wv ■ - &F sin A < 121 > 

where r 2 = s 2 + s, 2 - 2 s i; s i( cosy, sii 
V W v w 

come from the cosine and sine formulas, 
along the upper right wing panel is the 
itself and equation (121). Thus for t hi 
i s 



respectively. The total normal velocity 
sum of that induced by the wing on 
x-wing the induced angle in the wake 


a ( d> ) = — 
W v IT 


G'($i) d*, + 1 

COS^-COStJ) 7T 


(cosycos^T -n) G 1 ( 4> n ) d<j>i 
Icos^i-ncosy )2 + n 2 sin 2 y 


J o 


J 0 


where the dimensionless terms are 


( 122 ) 


cos* = n = s/s„, G = r/2s V (123) 

e g 

With the Fourier series for G in equation (1) and with equation (66), equation 
(122) simplifies to 


a - 
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z 

= 1 


sin n<j> 
s i n 4 > 


P ) 
n y 


n a 


(124) 


where 


P 

n 


1 

TT 


(cosy COS<h-r)) cos n<fr 7 d<f>i 

(cos^j-ncosy )2 + n 2 sin2 y 


J o 


(125) 


This integral is evaluated in equations (C2) 
P n integral is similar to the I integral of 
can be derived directly from equation (68) 
by a transformation of the variables. From 
the components of n v shown at the right, the 
relations for and are 


through (C8) in appendix C. The 
equation (69). Equation (124) 
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026) 


n v = m cosy = cosy cos^ , 5 = ni siny = siny cos^ 


These values of n y and 5 are substituted into equation (67) for cos^j and £, 
respectively, in the term which factors cos n^. Then I of equation (67) 
becomes P p as given in equation (125). This transformation method can be 
applied to any two or more spatial arrangements of separate planar vorticity 
sheets. 

The induced drag expression, derived in a similar manner to that of equa- 
tion (69), is 




n a 2 + 2 
n 




nn : 



where 


P 

nn 


* 


r 

1 

IT 

J 0 


P^* sin<f> sin n<j> d<j> 


027) 


(128) 


where P n * is given in equation (125), but with n = n*, and n = cos$. The drag 
influence coefficients, P nn *> are developed into a recurrence formula in 
appendix C, equations (C15) through ( Cl 9 ) . 

The lift coefficient given in equation (4) is with the lift normal to the 
wing from the aft view. For a cruciform wing or a tilted or rolled wing, the 
lift upwards is a component of the value given in equation (4), then 

Cj_ = 3-j cos J (129) 


where C L , A, and C Q are based on one wing of the crossed wings of the cruciform, 

i.e., the lift, induced drag, wing area, and width (-s g to $ e ), of one of the 
two equal wings. Thus 


( 2s e ) 2 L D. 

S 2 s e c av’ A S • C L qS’ C D i ’ qS 

For symmetric spanwise loading only odd integers n and n* 
equation (127) with (129) becomes 


(130) 

apply, then 
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oo a n oo a n * 
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(131) 


odd 
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a n /a-| coefficients for minimized induced drag. - Derivatives of equation 

(131) are obtained in the same procedure as that leading to equations (73) 
and (74). Then 

mP lm + P m! + ™ 57 + % ^ mP nm + " P mn } 57 = 0; m = 3 ’ 5 ’ 7 ’ • • • N < 132 > 

I n = 8 1 

odd 

where N is an odd integer. A simultaneous solution of these (N - l)/2 equations 
results in the evaluation of the a n /a^ ratios. These a n values, through 
equation (1), define the spanwise loading distribution for minimum induced 
drag for a cruciform wing where the wings are angled Y to each other in the 
yz plane. With the a^/ a ^ ratios known, then 1/e can be determined from equa- 
tion (131), induced angle in the wake from equation (124), wing bending moment 

ratio from equation (11), n from equation (12), and spanwise loading from 

cp 

equations (1 ) or (5) . 

an/a-j coefficients for minimized induced drag and wing bending moment, 
and antisymmetric loading. - The added constraint of minimized wing bending 
moment is taken into account identically as done in the derivation of equation 
(77). For the cruciform wing solution the I 1 s of equation (77) are replaced 
by P 's. Or the equations are simply, mtf m added to the left side of equation 

(132) m 

An antisymmetric loading solution can be obtained similarly as that lead- 
ing to the equation (81), except the I 1 s are replaced by P nm 's. With 
determined from equation (81), the drag efficiency factor is found from equa- 
tion (80) in which I is substituted by P nm - 
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Cruciform wing wi t h y = 90°. - In a cruciform wing when y = 90° there is 
no induced normal velocity at one vorticity sheet due to the other. The 
solution thus simplifies to the case of a single planar wing banked at an 
angle y/2. For the constraint condition of lift, the loading is elliptical 
along the span, n = s/s e . The C L is given in equation (129), then from equa- 
tion (131), e _1 = 2. For constraint conditions of lift and bending moment’ 
the loading characteristics are the same as those of equations (30) through 
(39), except that since by equation (129), C L (y) = cos (y/ 2), then is 
replaced by 2^ in equations (35) and (38), then 

C b = - kt), e" 1 = 2 (1 + kt 2 ) (133) 

When y = 90°, the influence coefficients P n and P^ m should be zero if 

there is no normal velocity induced from one vorticity sheet onto the other. 

Examination of equations (C2) and (Cl 5) shows that when y = 90°, the P n and 

P are zero for odd integer values of n and m, that is, for symmetric span- 
nm 

wise loading. For anti symmetrical loading, P n and P nm have values and thus 
the two wings are not aerodynamical ly independent of each other. 


Flow Field at a y Banked Plane, Induced by a Flat Vorticity Sheet 

The cruciform wing solution is also the solution for determining the 
flowed field about a banked plane, as 
shown at the right. At the polar coor- 
dinate point (n,y) the normal velocity 
to the y plane is given by the second -1 
term in equation (124). Thus 




"Vv 

(n,y) 


w(ri,y) _ 
V 


n P a 
n n 


(134) 


n = 1 

The lateral velocity along the y plane is analytically similar to equation 
(121), except that it is the cosine component, thus 

dv = cos A 
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where by the sine formula, cos A = (s v siny)/r. In dimensionless parameters, 
with equation (1), and integrating 


v_ 

V 




-sin y 


n 


2 

= 1 


n a — 
n it 


cosfticos n<h djj 
(cos^j-ncosy) 2 + n 2 sin 2 y 


(135) 


Since 2 cos^ cos n^ = cos (n + l)^ + cos (n - 1 ) <#> x , then equation (135) 
reduces to 


v(n»y) _ 
V 


sin y oo . 

- — p — ^ 2 n (J 

c n = 1 


n+1 


+ J ■, ) 
n-1 7 


(136) 


where the J n values are given in equations ( C 1 1 ) through (C13). The 
values of equation (134) are given in equations (C2) through (C8). Equations 
(134) and (136) are exact solution equations for predicting the normal and 
tangential induced velocities on a y-banked plane, induced by a flat vorticity 
sheet of arbitrary loading distribution, i.e., arbitrary a n (see eq. 1). 


V-Wing 


The V-wing induced angle equation is similar to that of the cruciform wing, 
except the spanwise integration for velocity normal to the surface is made 
only over the upper wing panels. Then equation (122) becomes 


a ( <J> ) = — 
W V TT 


r 7r /2 p t 

G'(*i) djjL + 1 

C0S<|)1- COScj) tt 


tt/2 


(cosycos^i-n) G 1 (jft] ) dft] 

(cos<|) 7 -ncosy) 2 +n 2 sin 2 y 


(137) 


By adding and subtracting the first integral with limits tt/2 to tt and noting 
that the pertinent geometric parameter for the V-wing is the dihedral angle 
related to y by 


r = ly 

1 V 

then equation (137) becomes 
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f 77 
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(cos ( {) 1 -n)[ (cos 4 >rncos 2r) 2 + n* sin 2 2rj 


r * 


tt/2 


(139) 
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With the G series of equation (1) 


<>n> n “n 


( 140 ) 


where 


nr 


= - — sin 2 r 


tt/2 


(cos^T+n) cos^t cos n<h 

(cos^-nHlcos^-nco^r) 2 + n 2 sin 2 2rJ 


(141) 


Using equation (140), the induced drag equation, derived similarly to 
that of equation (69), is 


V f ( n f 7 n V + 2 n ! 1 a n n J / V> 


(142) 


where 


fir 


Q nn * 

nn tt 


Q n * sin $ sin n$ d^ 


(143) 


where Q n * is given in equation (141) in which n is replaced by n*, and n = cos: 
The lift coefficient is 


C L = a-! cos r 


(144) 


where C^, A, and are based on S = 2s g c av , A = (2s e ) 2 /S, = L/qS, = 

D^./qS. For symmetric spanwise loading only odd integer n and n* apply, then 
equation (142) with (144) becomes 


ttAC 
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2 Z ^ S n* Q * 
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odd 
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Recurrence formula for Q n and Qnn*. - To formulate a recurrence formula 
for the induced velocity influence coefficient it is simpler to work with Q 
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in the form 


Q = q - r 
x n n 


(146) 


where 
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COS<f>! - 71 
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The recurrence formula for q and r are derived in the same manner as detailed 

^n n 

in appendices B and C. Then 
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Noting that q = q p and r = r , then 

q 2 = cos 2r - n + 4n cos 2r q 1 - (1 + 2n 2 ) q 
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q 3 = cos 2r + ^ n + 4n cos 2r q 2 - (3 + 4ti 2 ) q^ + 4n cos 2r q ( 
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The induced drag influence coefficient in terms of q and r, is equation 
(146) with n = n*, inserted into (143), then 


^nn* %n* " r nn 
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(153) 



where 
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Then using the recurrence formulas of equations (149) and (151) 
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a n /aT coefficients for minimized induced drag. - In equation (145), since 
a-j represents lift, this equation has the constraint condition of lift, 
that is, induced drag per unit lift. The partial derivative of equation (145) 
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with respect to a m /a-j is identically similar to that of equation (73) and (74). 
In equation (74) I mn is replaced by Q mn , then 

m ^lm + Q ml + m if + J 3 (mt) nm + ^ ‘ m - 3,5,7. . .N (160) 

odd 


A simultaneous solution of these equations leads to evaluation of the a n /a^ 
ratios. With the a^/a-j values known, then the spanwise loading is determined 
from equation (5). This is the loading distribution for minimum induced drag 
per unit lift of the V-wing. Other aerodynamic characteristics are 1/e from 
equation (145), induced angle normal to the wake from equation (140), wing 
bending moment ratio, cosr, from equation ( 11 ), and n C p from equation ( 12 ). 

The solution for the constraint conditions of lift plus bending moment is 
obtained from equation (160) which is modified by adding the term mtf , similar 
to that shown in equation (77). 


V-wing in ground effect. - When the apex of the V-wing touches ground the 
vorticity sheet system including ground images becomes that of the cruciform 
wing except the sense of the ground image vorticity is reversed. This can be 
taken into account by reversing the sign on the integrals in equation ( 122 ) 
when integrating along the lower panels of the sheet. Thus 
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With G 1 ^) = Ena cos n^, comparing with equation (124), the first 
n = 1 n 

integral of equation (161), has the term sin -n<j»/si n<#> , the second integral. 
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comparing with equation (125), has the term -P , and the third and fourth 
integrals, comparing with equations (146) through (148), has the term 2Q n . Then 
equation (161) becomes 


<* W U) = 


? (SinjT^ + 2 q - p ) n a 
v st n<j) y n n' n 


n = 1 


(162) 


The solution for the V-wing in ground effect follows the same procedure as 

that for the V-wing, but with the substitution of 2Q - P for Q , and 2Q * - 

n n n nn* 

P nn * for Q nn * in equations (141), (142), (145), and (160). The influence 
coefficients of P p and P nn *> in which y = 2r* are given in appendix C. 


Planar-Wing Winglet Configuration 

An aft view of this configuration consists of a flat wing span with wing- 
lets at the tips, as shown below. 



The dihedral cant angle, y, of the winglet is the angle between the extended 
wing arid the winglet. Theoretical aerodynamic characteristics of the planar- 
winglet have been developed in reference 9. The minimization method objective 
is to determine the loading distribution along the span extent of the wing and 
up the winglet such that induced drag for a given lift, and induced drag for 
a given root bending moment, are minimum. 

Planar-wing winglet aerody namic cha racteri st ics from reference 9. - From 
equation (69) of reference 9, the induced drag coefficient of a planar-wing 
winglet configuration with symmetric loading is 
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where (noting that ng = cos<})q) 
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For the winglet span position of <j)g = 5tt/32, or ^g = .881921 , values of 
D nn ^ had been computed for y = 90°, and 75°, and presented in table II of 
reference 9. This table of D nn * has been reproduced as table 13 in the present 
report. 

From equations (55) and (59) of reference (9), the lift coefficient and wing 
root bending moment coefficient are 
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Equations (163) and (169) are combined to form e _1 = tt ACg /C^ 2 , then 
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The ratio of equation (170) to (169) leads to 
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These induced drag, lift, and bending moment coefficients and aspect ratio 
are based on the extended span of the planar-wing span plus winglet spans, 
and on surface area of the planar-wing area plus the winglet areas. 

an/aT coef ficients for minimized induced drag fo r a given lift. - The 
partial derivative of equation (172) with respect to a m /a-| is 
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where m = 3, 5, 7, . . .N. The minimization solution is with equation (174) 
equal to zero. That is, the quantity within the braces in equation (174) must 
equal zero. Then 
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Equation (175) involves the solution of ( N- 1 ) / 2 nonlinear simultaneous equations. 
It is nonlinear because of e in the last term. A method of solution is to 
iterate e with estimates for e computed from equation (172), then equation 
(175) can be solved as linear simultaneous equations. With these new values 
of a /a.j , e is recomputed from equation (172), and the linear solution repeated. 

A first estimate for 1/e is that for elliptic loading, then a p /a^ = 0 for n > 1. 
Then from equation (172) for elliptic loading 
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a n /ai coefficients for minimized induced drag and wing root bending 
moment. - For this minimization process, equation (76) must be satisfied, thus 
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where the first term is given in equation (174). The partial derivative of 
equation (173) with respect to a m / a -j » inserted into equation (177) leads to 


2m — — + ~ D , + m -- D-, + E (m -- D + n -- D ) — + (1 - cosy) T 

a -j tt ml 7T 1m n= ^ t nm tt mn a-j -ire m 


odd 
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m+1 


[1 - I 0 - ^sr) (T-| - T n ^)] + 8t } 0 - cosy) T, 


n=3 

odd 
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cos<f) 0 ] [1 - - (1 - cosy) (T ] 
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odd 



- cosy) T n 


cos4>q] 



1 

2 


0 


m = 3, 5, 7, . . . N 


(178) 


where t is a constant similar to that developed in equation (20). The solution 
of equation (178) with respect to a n / a -j is the same as that of equation (175), 
but with the addition of the t term to the linear solution as e is iterated. 

Example numerical solution with = 5 tt/ 32. - With the D nn * for y = 9CP , 
and T , of tables 13 and 14, inserted into equation (175), this equation simpli- 
fies to the numerical form, for m = 3, 5, and 7 
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.001665 + 2.821910 - .166305 - .016025 — + .117625 ~ = 0 

S'! 3 ^ Q -j G 

-.002125 - .166305 — + 4.908620 P + .099925 P + .126340 f- = 0 

a-j a-j a-j e 

.011935 - .016025 — + .099925 — + 7.184340 ~ + .076970 - = 0 

a-j a-j a-j e 

where 


(179) 


3 o 3 r 3 -j 

E = .952165 - .117625 — - .126340 - .076970 -7- 

a l a l a l 


(180) 


Equation (172) becomes 


1 = tL [1.03335 + .00333 — - .00425 — + .02387 - +2.82191 (— ) 2 + 

e E z L a-j a-j a-j v a^ 

3 3 3 3 3 3 

4.90862 (— ) 2 + 7.18434 (— ) 2 - .33261 — — - .03205 — — + 

a-| 3-| 3-j 3-j 3-j a-j 

3 3 

.19985 (181) 

a l a l 


and equation (173) becomes 

C. = nr (.30020 + .11853 ^ - .13513 - .03109 -J-) (182) 

D Trt 3 -j 3-| 3 -j 

For elliptic loading, equations (176) or (181), and (182) become 


— = 1.139785, C. = .100357 (183) 

e be 

c 

A one equation solution is obtained from the first equation of equation (179), 
starting with the e value of equation (183) and iterating until e converges. 
Similarly for a two equation solution, but starting with the 1 - eq. solution 
of e, etc. The converged three through one equation solutions are listed as 
follows : 
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1 

a 3 

a 5 

a 7 

c 


e 

a l 

a l 

a l 

c b 

3 eq. solution 

1.127466 

-.047592 

-.028843 

-.012991 

.098846 

2 eq. solution 

1 .128776 

-.047542 

-.029110 


.098826 

1 eq. solution 

1.133285 

-.045823 



.097987 


The three equation result is an accurate solution for the wing with a winglet 
at <j>Q = 5tt/ 32, or hq " .881921 , with y = 90 degrees. The loading is given by 
inserting these a n /a-| values into equation (85). This loading is along the 
span extent of the wing and up the winglet to the winglet tip. The spanwise 
loading for minimum induced drag for a given lift is inboard from elliptic 
loading, particularly the loading is substantially less over the winglet as 
indicated by the negative values of a n /a.j . Compared with elliptic loading 
values the minimized solution has a -1.081 percent smaller e” 1 , and a -1.506 
percent smaller C^. The e _1 value of 1.127466 is with C^, , and A based on 

the combined areas and spans of the planar wing and winglet. When based only 
on the planar wing area and span, this e " 1 value is multiplied by hq 2 > then 
e " 1 = .876926. This means that the addition of two 13.39 percent wing semi- 
span winglets at y = 90 degrees, to a planar wing, can reduce e " 1 from one or 
larger to .877. 

The equations for the minimized induced drag and bending moment solution 
are given in equation (178). These equations are the same as those of equation 
(175) but with an addition of a t term. For 4 ^ = 5n/32, y = 90 degrees, the t 
terms to add to equation (179) are 


t (.592669 - .123477 ■— - .051120 ~) 

a l a l 

3o 3 7 

t (-.362955 + .123477 + .025891 — ) > 

a i a i 

Si o 9 p 

t (-.025992 + .051120 — - .025891 — ) 

a l a l J 

For t = 1 in equation (185) and adding to equation (179), gives 


(185) 
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a 3 a 5 a 7 E 

.594334 + 2.82191 — - .289782 t 2 - - .067145 + .117625 ^ = 0 

d-j d-j d-j 6 

-.365080 - .042828 — + 4.90862 — + .125816 ~ + .126340 ~ = 0 

cl-j a-j a-j 6 

a 3 a 5 a 7 F 

-.014057 + .035095 — + .074034 — + 7.18434 ~ + .076970 |- = 0 

a i a i a i e 


> 0 86 ) 


where E and 1/e are given in equations (180) and (181). In the iteration of e 
it is simpler to do a one equation solution of equation (186) for the first 
estimate of e. The three equation iterative solution of equation (186) gives 


3^ aj, di-j 

— = -.25932, — = .03994, — = -.01069 

a l a l a l 


- = 1.2887, C. = .086015 
e b 


087) 


where is determined from equation (182). The spanwise loading is given in 
equation (85) which with these a n /a -j coefficients shows that this loading is 
strongly inboard of elliptic loading. 

Combining equations (183) and (187), the e and ratios are 


— = 1 .13065, 
e 


'be 


= .85709 


(188) 


From equations (49) and (50), for the same lift and wing root bending moment 


s C, 

— = ^ = 1.16674 (189) 

S Li 

ec b 

D. 2 e„ 

!=(_§£) — = .83058 (190) 

D. s rt e 

ic e 

where subscript c denotes a configuration with elliptic loading, and s e is the 
semispan of the planar wing plus the span length or height of the winglet. 
Equations (189) and (190) means that for <j>g = 5tt/32, y = 90 degrees, if the 
wing with the inboard loading of equation (187), has a span length 16.67 per- 
cent longer than the wing of elliptic loading, then the induced drag will be 
17 percent less, for the same lift and the same wing root bending moment. 
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RESULTS AND DISCUSSION 


In the present theory, minimization solutions are obtained by expressing 
the spanwise loading in a Fourier series, determining the induced angle normal 
to the wake, formulating the induced drag as a series of products of Fourier 
loading coefficients and induced drag influence coefficients, imposing con- 
straint conditions, then minimizing the constraint conditioned induced drag with 
respect to the Fourier loading coefficients. This results in a set of unknown 
Fourier loading coefficients with an equal number of equations. The solution 
of these simultaneous equations leads to values of the Fourier loading coef- 
ficients and thus the loading for minimum induced drag. 


Planar Wings 

Constraint that either lift or rolling moment is specified. - Example 
solutions are shown at the beginning of the planar wing chapter. For the 
constraint that lift is specified the solution is 


°i C _ 4 . , „ _ C L 2 „ _ f l _ 4 Mn-M 

C L C av 77 S1P C D i irA 9 6 1 ’ C b 4 77 * \p~ 3 77 ^ 19 ^ 

where sin<£ = (1 - ti 2 ) 2 , i.e., elliptic, and f-j is given in equations (9), 

(10a), or (36). 

For the constraint condition that rolling moment is specified, the solution 
for a planar wing i s 


c o c 1 6 

C D 


C £ c av 


32C £ 2 

tt/T" * e r = 1# n cp 


3 77 

16 


(192) 


no- 
where sin 2<j) = 2 n( 1 - n 2 ) 2 , and C £ is wing rolling moment coefficient. 

Constrain ts of lift and of wing bending moment about nb« - The minimization 
solution in terms of infinite series is presented in equations (21) through 
(26). The solution is in terms of a constant, t, which governs the quantity 
of the bending moment constraint. When t is zero the solution simplifies to 
that of equation (191). An upper limit for t is given in equation (29) at 
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which point the loading becomes negative at the wing tips when t > t-j . Series 
summation methods are developed in appendix A so that equations (21) through 
(26) can also be represented in terms of closed functions. In terms of closed 
functions the planar wing solution of minimized induced drag for a given lift 
and a given wing bending moment about span station is presented in equations 
(30) through (39). For the constraint condition of wing root bending moment, 
tk = 0, then equations (30) through (38) simplify to equations (40) through 
(44). For n b = 0, equation (41) can be written as t = 6 - 9-n-n C p/2, which 
substituted into equations (40) and (44), results in a loading function and 
induced drag function which correlate exactly with those presented in reference 
2 which is a wing-root bending moment method. Equations (30) through (39) 
represent an independent solution as determined by the use of the present 
minimization theory, and developed originally without knowledge of the previous 
work of references 1 through 4. In terms of n^, equations (30) through (39) 
are unique in the literature. 

Numerical values of k, k-j , f -| , and kg for various values of are 
presented in tables 1 and 2. With these constants the loading distribution, 
spanwise center of pressure, induced angle in the wake, wing bending moment 
parameter, and induced drag efficiency factor, are given in equations (30), 

(32), (34), (35) and (38), respectively, for arbitrary values of t. From 
equation (30) it can be seen that the loading distribution is elliptic loading 
times the factor, 1 - kgt. With positive values of t, and since in table 2 
kg is negative inboard and positive outboard on the span, then this optimum 
loading will be inboard of elliptic loading. The aerodynamic characteristics 
for specified values of t are given in tables 4 through 8. Greater detail is 
given in the section 'Numerical values of loading characteristics for various 
n b and t‘. A measure of the efficiency of wings with inboard loading is to 
compare induced drag with a wing with elliptic loading having the same lift 
and same wing bending moment. This is done in the section entitled 'Comparison 
with Elliptic Loading'. Aerodynamic characteristics ratioed to those of the 
wing with elliptic loading, for same are presented in table 9. Correspond- 
ing spanwise loadings are given in table 8. In addition the data of tables 4 
through 7 contain some comparisons, that is, the C^ 2 / e ratios are the same as 
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V D ic* and C bc /C b = b / b c- 

Tables 10 and 11 contain data for the condition that wing bending moments 
are equal about the same y b = y bc span station. Since the elliptic loading 
span is smaller, then when y^ = y bc , the dimensionless coordinates become, 
n b > n bc* sma ^ er % favors a higher bending moment which means there 

must be a larger inboard shift of loading which, in turn, increases induced drag. 
The induced drag ratios of table 10 are higher than those of table 9 where 
n b = n bc‘ Corresponding spanwise loading coefficients are presented in table 
11. When b/b c = 1, the loading is elliptic, that is 4(1 -d 2 )^/^ The loading 
for the asterisk marked lines in table 10 is wingtip zero-slope loading pre- 
sented at the top of table 8. 

Derivation of t for minimal induced drag and wing span. - Shown in the 
data of table 9 is a relative decrease of induced drag with increase in wing 
span, for wings with constant lift and constant wing bending moment. A design 
objective could be to get the smallest induced drag for the least increase in 
span. That is, to minimize the product of the span ratio and the induced drag 
ratio. The product of equations (52) and (57) gives, for same n b 


bD i 

b c D ic 


* 


= (1 + kt 2 ) (l - 4- t) 
T 1 


(193) 


where e is from equation (38). Equation (193) is minimized by taking the 
partial derivative with respect to t and equating to zero. This produces 


t = 3^ [1 - (1 - 


3k 

TJ 


1 


(194) 


b_ = 1 = 3 

b c 1 - ~ t 2 + (1 - ~ T )% 
f l f l 

1 = 1 + kt2 = | + [1 - (1 - fV 8 ] 


(195) 

(196) 
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(197) 


?r— = l (£-)' 2 = (1 + ktz) (1 - f- t) 2 
u ic e D c T 1 


Equations (194) through (197) are valid when the span-drag product of equation 
(193) is minimum, and when both wings have the same n^. 

For the constraint condition of wing-root bending moment, = 0. Then 

from equations (36) and (37), or table 1, k = 2/9, and f-j = 4/3. Using 
equations (193) through (197) gives 


t - .418861 



1.075049 


e = 


C L 2 


ttAC 


D i 


.962475 


= .898987, 
u ic 


b D. 


D ic 


.966456 


(198) 


Using the above t, the spanwise loading distribution, and nc p> a w , and C^, are 
determined from equations (40), (41), (42), and (43), respectively . The values 
of equation (198) show that for a wing with the loading of equation (40), with 
just a 7.5 percent increase in span, a 10.1 percent decrease in induced drag 
can be realized compared to a wing with elliptic loading, with both wings at 
same lift and same wing root bending moment. 


For the constraint condition of wing bending moment about span station, 

= .2, then k = .176851 , f-| = .784747, which are obtained from equations (36) 
and (37), or from table 1. Using equations (193) through (197) gives 


t = .928705 


e = 


C 2 


irAC 


D i 


.867654 


> (199) 


= 1 .264693 


D. b D. 

= .720582, t-— 4 — = .911315 

u ic D c u ic 


Using the above t and the spanwise loading distribution, and n C p, a w > and 
C^, are determined from equations (30), (32), (34), and (35), respectively . The 
values of equation (199) show that for a wing with the loading of equation (30), 
with a 26.47 percent increase in span, a 27.94 percent decrease in induced drag 
can be realized compared to a wing with elliptic loading, with both wings at 
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the same lift and wing bending moment about span station n b - .2 = n bc * Data 
from table 9 for n b = .2 show that with just a 7.5 percent increase in span, 
a 11.8 percent decrease in induced drag can be realized. 

Nonplanar Wings 

Example applications of the minimization theory are made to the config- 
urations described as, biplane, wing in ground effect, cruciform wing, V-wing, 
V-wing in ground effect, and planar-wing winglet. The biplane solution is 
obtained from equation (74). This equation is valid for a biplane with wings 
of equal span and same spanwise loading distribution; is independent of stagger; 
and is only dependent on height between wings or gap since I nn * is a function 
of c. Equation (74) is also valid for the case of a planar wing in ground 
effect, provided that I is substituted by -I nm , and that C = 4h/b, where h 

is height from ground. For s = 1/2, the solution for minimum induced drag for 

a given lift, for the biplane, the loading characteristics are those given in 
equations (84) and (88), which show an outboard shift of loading compared with 
elliptic. Similarly, for s = 1/2, or h/b = 1/8, for the wing in ground effect, 
the loading characteristics are those given in equations (97) and (98), which 
show an inboard shift of loading compared with elliptic. For C = 1/2, the 
solution for minimum induced drag for a given lift and a given wing root bending 
moment, for the biplane, the loading characteristics are given in equations (90) 
through (96), while for the wing in ground effect, in equations (99) through 

(103). Compared to a wing with elliptic loading, the biplane with a 14.1 per- 

cent greater span, has 12.8 percent less induced drag, while the wing in ground 
effect with 28.6 percent greater span, has 13.6 percent less induced drag. 

The minimization theory application to a cruciform wing is similar to 

that for the biplane. The primary difference is that the induced drag influence 

coefficient, P ., replaces I * of the biplane. Minimized induced drag solu- 
nn* nn* 

tions are obtained from either equation (132) or equation (132) plus mtf m , which 
are analogous to equations (74) and (77), respecti vely . The P nn * integrals are 
evaluated in appendix C. 

The analyses for a V-wing is the same as that for the cruciform wing except 
that the spanwise integrations for the induced velocity influence coefficient 
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is for half the wing span. The induced drag influence coefficient is denoted 
by Q nn *- Minimized induced drag solutions are obtained from either equation 
(160) or equation (160) plus mtf^, which equations are analogous to equations 
(74) and (77) with I nn *> or equation (132) with P nn *- When the pointed end of 
the V-wing touches ground, the vortex sheet plus ground image vorticity resem- 
bles in part the cruciform wing (see eq. 161). The solution for the V-wing in 
ground effect (V-wing apex touching ground) follows the same procedure as that 
for the V-wing, but with the substitution of 2C) - P for Q. and 2Q - P 

for Q nn * in equations (141), (142), (145), and (160). The P n and P nn * coef- 
ficients, in which y = 2r, are given in appendix C. 

The planar-wing winglet configuration influence coefficients are derived 
in reference 9 and reproduced in the present report in equations (164) through 
(168) and in table 13. The minimization solution is obtained from equations 
(175) for a given lift, and (178) for a given lift and a given wing-root bend- 
ing moment. The numerical solution for <j>q = 5tt/32 (winglets starting at n g = 
±.88192) with winglets at right angle to the wing, is presented in equation 
(184) for the constraint of a given lift. This optimized loading is inboard 
relative to elliptic loading, particularly the loading is less than elliptical 
over the winglet. Comparing the same planar wing with and without added wing- 
lets, the induced drag is 12.3 percent less for the wing with winglets. The 
solution with constraints of lift and wing-root bending moment is presented in 
equation (187). The results show that if a wing-winglet with the optimized 
inboard loading has a span length 16.7 percent longer than that of the wing- 
winglet with elliptic loading, then the induced drag will be 17 percent less, 
for the same lift and wing-root bending moment. 

Other Applications of the Analysis 

An exact solution for the flow field about a flat vorticity sheet or 
planar wake is presented in equations (118) and (119). These equations are 
the vertical velocity or downwash, and lateral velocity or sidewash, induced 
by a flat vorticity sheet of arbitrary spanwise vorticity or loading distri- 
bution. The loading distribution is represented in equation (1) in which the 
a n are Fourier loading coefficients which can have arbitrary values. Similarly, 
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equations (134) and (136) are presented for predicting the normal and tangential 
induced velocities on a y-banked plane, induced by a flat vorticity sheet of 
arbitrary loading distribution. 

With equations (118) and (119), the induced velocities can be predicted 
at any point in space. These equations can thus be used to set up the mini- 
mization solution for such configurations as multi -wings, wings in formation 
flying in arbitrary pattern, and banked wings. The minimization equation for 
the V-wing in ground effect is given in equation (162). The V-wing with apex 
at the boundary of an open wind tunnel has an identical solution to the 
cruciform wing presented in the section Cruciform Wing. The biplane solution 
is also the solution for wing flying near a wake, that is, over or under a wake 
which has the same loading distribution as that of the wing. The minimum induced 
drag solution for a wing under an arbitrarily loaded wake, is obtained by fol- 
lowing the procedure in the biplane section. Let the wake loading be 


G 


w 



b n sin n<j> 


( 200 ) 


then 


°w 


£ n a + 
n = -j sin <p n 


I n b 
n n 


n = 1 

With equation (201), 1/e of the wing becomes 

C 

V -) 2 + ? 
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( 201 ) 


n 


? n (V + 2 . 

n = 1 a l l L n = 1 a l n*= 1 


n* I 


nn* b 


( 202 ) 


The a n /a -j coefficients for minimized induced drag are obtained by taking the 
derivative of 1/e with respect to a m /a-j, and setting equal to zero. This 
solution gives 


a 


n 


a 


1 



oo 

£ 


n* = 1 




(203) 


where the lift coefficient ratio is that of the wake loading to that of the 
wing. The values of equation (203) in equation (1) give the loading distribution 
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of the wing with minimum induced drag, flying at a ^ distance above or below 
a wake which has the loading of equation (200). When the wake loading distribu- 
tion is elliptic, then b p = 0 for n > 1, and equation (203) reduces to 


a 


n 


a 


1 




n 


(204) 


which loading is outboard relative to elliptic loading. 

The following two sections contain other applications in greater detail. 


Flow Field Solution of a Thin Wing Chordwise Vortex Sheet 


Two-dimensional thin wing theory and wake theory are in principle the 
same, differing only by coordinate definitions. In equations (114) and (115), 
substitute the longitudinal coordinate x for y, y for dr/dy-j , c for b, and 
~u s ( x , z ) for v s (y,z). Then 


w (x , z ) = — - 

S Ztt 

u_(x,z) = I- 
S 


'c (x-x-j ) y dx-| 

. 0 (x-x 1 )2 + Z2 

fc Y dx l 
J 0 (x-X-,)2 + Z2 


(205) 

(206) 


where y is circulation per unit length along the wing chord, related to total 
circulation and local lift coefficient by 




J o 


(207) 


In terms of dimensionless coordinates, chordwise trigonometric coordinates, and 
chordwise loading represented by a cotangent term plus Fourier series, define 


| = ? = j (1 - cose), | = S 2 (208) 

= - \ AC p = 2 A q cot Jr + 2 z A n sin ne (209) 


With equations (208) and (209), equation (205) becomes 
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w s (£,fi) -j 


[cos^->(1-2^)]|a o (1+cos0-j )+j e A n [cos(n-l )e-j -cos(n+l )e-| ]|d 0 i 

[cose-] -~(T - 2c )T 2 + 4’si 2 


( 210 ) 


These integrals compare with the integrals of equations (67) or (Bl), then 


w s (c.fl) = A 
V M ( 


< I o + V 


*1 


E 
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A„ (I 


n-1 


- W 


( 211 ) 


where I = I (1 - 2c, 2n). By a similar derivation, noting correspondence 
with the integral of equation (B12), equation (206) becomes 


u s (c»fi) 


= 2 [A„ (J, + 


J l ) + 2 


E 

= 1 


A. (J 


n-1 


^n+l 


( 212 ) 


where J n = (1 - 2c , 2n) . These I n and J n induced velocity influence coef- 

ficients are evaluated by substituting 


ti = 1 - 2c, and ^ = 2ft (213) 

into equations (B7) through (Bll), and into equations (B13), (B14) and (B18) 

through (B20). With these n and substitutions the p and r values of I and 

n 

J n becomes (from eq. Bll) 

P = 4(c 2 + ft 2 ), r = 4 [ ( 1 - c) 2 + ^ 2 ] (214) 

The dimensionless coordinate c is measured from the leading edge of the chord- 
wise vortex sheet, and ft is measured from the sheet. 


Equations (211) and (212) represent an exact analytical solution for 
induced downwash and induced longitudinal velocity, due to arbitrary chordwise 
loading as represented in equation (209). The total longitudinal velocity is 



u s (c>a) 

V 


(215) 


where V is freestream velocity. 
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Formation Flying of Wings with Wingtips Linked 

The objective is to link wings together and determine the spanwise loading 
such that induced drag is minimum with the constraint conditions that the bend- 
ing moments at the wing connections be zero. This is because the term linked 
here means hinged and cannot support a 
bending moment. A three wing symmetri- 

i i — i . i 

cal combination shown here, has wings 
spanning the distances -1 to -n b , ~n b 


-1 


-n, 


0 


1 


to n b , and nb to 1 . 


The constraint c Pb 

condition when no lift transfers through the link is that the center of pressure 
of the outer wing loading, n C p 5 » be at the midpoint of the nb to 1 span. 

The minimization of induced drag for a given lift and wing bending moment 
has already been developed in the present planar wing chapter. These results are 
valid for an arbitrary value of t. For this linked wing problem, t must be 
determined such that the constraint condition on n C p b is satisfied. The span 
position, Ti C p b > is the sum of n b and the ratio of bending moment about n b » to 
wing lift from n b to 1 . Then 

2C. 

ilk + 1 7T (216) 


'cpb 


L tA 


where L b /L is the fraction of lift of the outer wing compared to the total lift 
of all the wings. Using the relation a n /a-j = -tf n from equation (20), this 
ratio is 


where 


Gdrt 


Gdi 


~i 


h i = 




- t £ 
n = 3 
odd 


h f ) 
n n 


(217) 


4> b - sin <j> b cos <j> b , <j> b = cos 


(218) 


h = sin (n-1)i>h _ sin (n+l)4>b 
n n n-1 n+1 


(219) 
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The constraint condition is that n C p b = (1 + n b )/2* 
equation (24) for C b , and equation (217), leads to 


f T - 0 - V 

k - (1 - n b ) 


'1 


v h f 
L ~ n n 
i = 3 

odd 


Then equation (216), with 


( 220 ) 


where f -j , f n , and k are in equations (9), (10), and (37), respectively, or in 
tables 1 and 3. 


Equation (220) represents the minimized induced drag solution for three 
symmetrically linked wings with arbitrary outer wing span. With t determined 
for a given n b * the three-wing aerodynamic characteristics are evaluated from 
equations (30) through (38), and the lift ratio from equation (217). In the 
three-wing problem when lift is able or allowed to transfer through the links, 
then n C p b can be reduced until the net three-wing optimum loading is elliptical 
and net e is unity. 


For the case of n b = 0, the solution reduces to a linked two-wing config- 
uration. When n b = 0, then k = 2/9, k-| = 1/6, f-j = 4/3, h-| = tt/ 2, h p = 0 for 
n > 1, and equation (220) becomes 


t = - | (3 tt - 8) = -1.068583 (221 ) 

then the two-wing aerodynamic characteristics from equations (40) through (44), 
and (217) become 


■TT~ = I - *)(1 - n 2 ) % + (3ir 
C L c av " 

- 1 
n cp 2 

f£-a w = - (3 tt-9) + ( 3tt- 8) | n | 

= -.42478 + 3.3 57 05 [ n | 


- 8)n 2 


cosh" 1 




e = [1 + ^ ( 3 tt -8 ) 2 ] 1 = 

^b _ 1_ 

L “ 2 


> ( 222 ) 


.79761 


Thus, two equal wings flying linked together with the optimal net spanwise 
loading in equation (222), with no rolling moments will have 0.6269 as much 
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induced drag as when the wings are flying separately with elliptic loading. 

CONCLUSIONS 

An analytical minimization theory has been developed which is applicable 
to complex aircraft nonplanar configurations, with simple or complex constraint 
conditions. The application of the theory to many configurations and flow 
conditions have proven to be not excessively complicated and show that the 
analyses remain mathematically manageable and viable. This is because the 
solution method is based on a Fourier series representation which, in effect, 
breaks up ultra complex expressions into a sequence of analytically workable 
terms. The solution remains analytically exact since a Fourier series can be 
carried out to convergence, or final terms can be summed. In the planar wing 
minimization solution, the spanwise loading distribution is determined for which 
the induced drag is minimum for a given lift and wing bending moment about a 
given span station. The planar wing solution aerodynamic characteristics are 
given in equations (30) through (38). These equations are unique. They also 
provide the solution for three wingti p-1 inked wings. Minimization solutions 
are made for nonplanar type configurations. These show the optimum spanwise 
loading for different constraint conditions. In many cases the reduction of 
induced drag can be substantial, particularly with either the constraint of 
wing bending moment, ground effect, or formation flying. Some of the solutions 
show marginal reductions in induced drag, relative to that due to elliptic 
loading, but do provide the analytical loading for the minimum. 

For two of the simplified cases or conditions, previous theoretical sol- 
utions had been obtained by different analyses, which compare exactly with the 
present results. These include the planar wing with constraints of lift and 
wing-root bending moment (r^ = 0), and the flow field of a planar vorticity 
sheet with elliptic loading distribution. This minimization theory has proved 
to be a useful analytical tool for obtaining the exemplary minimization solutions 
in the present report. The method can be applied to many other problems of 
interest, such as, complex aircraft configurations, wings with arbitrary number 
of structural or performance constraint conditions, and optimized pattern for 
wings in flying formation. 
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APPENDIX A 


MATHEMATICS OF TRIGONOMETRIC SUMMATIONS 
Algebraic Series 

Need to evaluate series summations of series with terms of the type 


Z 


n 


n (n + a) (n + b) . . . (n + m) 


(Al) 


where n are odd numbered integers and m are arbitrary numbers. This family 
of infinite series begins with the series of tanh _ 1 Z S thus 

1+7 73 75 “ 7 *^ 

tanh l l = h In = Z + . . . = Z + n ~ 

1 - Z 35 n - 3 n 

odd 


then 


00 7 n+a 

z L 

n = 3 n 

odd 


-Z a+1 + Z a tanh _1 Z 


(A2) 


Similarly 


tan ~ l l = Z 


il+ ll 

3 5 


Z + £ , 

n = 3 

odd 


n -1 

(-1) 2 Z n 


then 

n -1 

j (-1) 2 Z n+a 
n = 3 n 
odd 



+ Z a tan _1 Z 


(A3) 


The l/n(n+a) terms summation is obtained by dividing equations (A2) and (A3) by 
Z, then integrating from 0 to Z, thus 
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°° -,n+a+b 7 a+b+l 7 a+b 7 b f 

■ c- W + V tanh ' lz -f- 


n = 3 n(n+a ) 
odd 


a ? 0 / - (odd integers) 


1 - Z2 


dZ] , 


°° 7 n+b 

z « ‘=— 7 ~ 
n = 3 n z 

odd 


[-z b+1 + z b 


tanh _1 Z 


dZ] , a 
0 


= 0 


(A4) 

(A5) 


00 7 n+b l . ■) 7 b 7 b 

n S 3 nfri^TT = Z - §- (1 + Z) In (1 + Z) + £- (1 - Z) In (1 - Z) 
odd 

a = -1 


n-1 


°° / t \ 2 7 n+a+b 7 a+b+l 7 a+b 

e l = r- - + - 

n = 3 n(n+a) L a + 1 a 

odd 


- = [- + 
a f 0 f - (odd integers) 


tan -1 Z - f- 

a 


1 +Z2 


0 


(A6) 


r Z a dZ ^ 


(A7) 


n-1 

l (-n 2 Z n+b 
n = 3 r)2 

odd 


[-Z b+1 + Z b 


tan ~ l l 


dZ] , a = 0 
0 


(AS) 


n-1 

? H) 2 z" +b 


n - 3 
odd 


n(n-l) 


= Z b+1 - Z b tan ~ 1 Z - In (1 + Z 2 ) 


a = -1 


(A9) 


Apply equations (A4) and (A7) for various values of a, leads to 
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(1 + Z) In (1 + Z) + 


oo 

E 

^n+b 

Z b 

n = 3 n(n+l ) 
odd 

2 

0 - 

z) In (1 

- Z) 

co 

Z 

^n+b 

y 1 
_Q 
M 

1 

n = 3 
odd 

n(n+ 2 ) 

2 

oo 

y. 

z n+b 

z b + i 

n = 3 n(n-2) 
odd 

2 

CO 

z 

n -1 

(-D 2 

^n+b 

n = 3 
odd 

n(n+l ) 

CO 

y 

n-1 

(-D 2 

z n+b 

n = 3 
odd 

n(n+2) 

oo 

E 

n-1 

H ) 2 

z n+b 


,b-l 


yU T I 7 u-z. 

^ 3 — - (1 - Z 2 ) tanh _1 Z 


7 b+l . 7 b-l 

+ Z tan _1 Z - ~ — ln 0 + Z 2 ) 


7 b-l 7 b+l 7 b-2 

~Z ~3 — + ~2 — ^ + Z 2 ) tan 


z b+i z b 


n = 3 
odd 


(A10) 

(All) 

(A1 2) 

(A13) 

(A14) 

(A15) 


In equations (A 6 ), (A9), (A10), and (A13), let b * 0 then taking derivatives 
with respect to Z leads to 

oo 7 n+b 7 b+l 

n= Z 3ir = -V ln 0 - Z 2 ) (A16) 

odd 

n = 3 n+T ' -T-- V 1 " 0 - Z2 ) (A17) 

odd 
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oo 



odd 


n-1 


Ml 


2 z n+b 
n-1 


n-1 

s (-1 ) 2 Z n+b 
n = 3 n+1 

odd 


7 b+l 

^2 — In (1 + Z 2 ) 

7 b+l 7 b-l 

V” + V” ln (1 + z2) 


Extension of equation (A2) leads to the two summations 


00 7 n+b 

z 

n = 3 n+2 

odd 



+ Z b " 2 tanh _1 Z 


00 7 n+b ,, 0 

H r> ~ ~o = Z b tanh ~ 1 Z 
n - 3 n-2 

odd 


The integrations of equations (A4), (A6), (A7), and (A9) result in 
summations where a and b have values such that no denominator 


“ 7 n+a+b+c 7 a+b+c+l 

2 ^ - r ^ i 

n = 3 n(n+aT(n + a+b) L (a+1 ) (a+b+1 ) 
odd 


7 a+b+c 

a[a+b]~ tanh " 1Z - 


Z c 

ala+bj 


Z a+b dZ 
1-Z* 



([ dZ)dZ] 

J 0 


E Z n+b+C = r b+2 7 b+c+l 

n = 3 n(n-l ) (n+b) L ( b+1 ) 2 L 
odd 


In (1 


Z 2 ) + 


Z c 

b(b+l) 


f b Z 

Z^ 

J 




yb+C 

V — tanh ~ l l 


^b+c+l 

2TB+TT 


(A18) 


(A19) 


(A20) 


(A21 ) 

the following 
vanishes: 


(A22) 


(A23) 
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« 7 n+c 7 c 7 ( 

n = 3 n(n-l)2 = ^ 2~ C + z ) ln (1 + z ) - J 
odd 


(1 - Z) ln (1 - Z) - 


,c+l 


IlL(Wfl dZ 5 

L- r\ 


(A24) 


n-1 

03 f n \ 2 ^n+a+b+c 

z 1 i L 


n = 3 n(n+a) (n+a+b) 
odd 


,a+b+c+l 


a+b+c 


-jKX ' U' L ' I ^ U « L/ ' O ^ 

= C ‘ (a+l )('a+b+Ty + FfS+bT tan ’* Z ' 


(a+bT 


7 a+b y c 
T+Z2" dZ ~ a 


,b-l 


' Z a dZ 
1+Z 2 


) dZ] 


(A25) 


n-1 

°° / x 2 7 n+b+c 

z v~ 1 \ . 4 

n = 3 n(n-lT(n+b) 
odd 


= c 


In (1 + Z 2 ) + 


b(b+TT 


b+2 -jb+c+l Z 
(b+1 )2 L 

b Z 
dZ] 

0 


1+Z2 


b+c 


,b+c+l 


b tan 1 1 “ 2fb+l ] 


(A26) 


n-1 

” (-1) 2 z n+b 

n = 3 n(n-l )2 
odd 


[Z b+1 ln Z - Z b+1 


tan _1 Z Z 

— ^ d Z — 


b+1 


l n li t* 2 . ) - dz j 
L 0 


(A27) 


Series of the type given in equations (A22) through (A27) can be determined 
by using these equations or by solution with partial fractions. Thus to 
evaluate the example 
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00 2 n+ 1 
n = 3 n(n-2) (n+1 ) 
odd 


then in equation (A22), a = -2, b = 3, c = 0, hence the summation is 


[ | tanh _1 Z + h 


z_dz + ^ 

1-Z2 * 


z 2 ( 


dZ 


Z 2 (l-Z z / dZ ^ 


thus 


“ Z n+1 Z; 
n = 3 n(n-2) (n+1 ) 3 
odd 


(1 - \ - §^) In (1 - Z) - 


-Mr 1 (1 + t - I-) In (1 + Z) 


(A28) 


By partial fractions, 

1 = 1 1 
n(n-2) (n+1 ) “ 3n(n-2) ” 3n(n+l ) 

then the summation is 1/3 of the difference between equations (A12) and (A10), 
and equals that given in equation (A28). 

Trigonometric Series 

The trigonometric series can be obtained from the algebraic series by 
letting Z be a complex function. That is 

Z = e 1 ^ = cos <j> + i sin <f> 

_,n in<J> . 

Z = e = cos nj>+ i sin n<f> 
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thus the cosine series is the real part and the sine series the imaginary part 
of the algebraic series. The real and the imaginary part of the functions in 
the algebraic series are listed as follows where absolute values are taken in 
the log terms: 


In (1 + Z) = In 2 cos | + i | 

In (1 - Z) = In 2 sin | + i (- + f) 

In (1 + Z 2 ) = In 2 coscf) + i $ 


In (1 - Z 2 ) = In 2 sincf) + i (- A + ) 

L. (p 


tanh ~ 1 Z = % tanh _1 cosep + i 

H ( }> 


tan ~ 1 Z = + i h tanh _1 sin<j> 

4 ( p 


In addition since dZ/Z = id<f>, then with equation (A30) 


(A30) 


tanh _1 Z 


( 2 2 " 1 - 1 - 1 ) B m 


dZ - - 5* + i[ f (in 2 + 1 - in*) - % m \ , ^ ^ jj). 


2m+l 


(A31) 


where B m are Bernoulli's numbers, 1/6, 1/30, 1/42, 1/30, 5/66, 691/2730, 7/6, 


tan ~ 1 Z 


(2 2m-1 - 1) B m 


dZ “ ^ - ♦) On 2 - 1 - in (J - *)] - * m S , 

(A32) 


(£- *) 2m+1 + 1 = 


72 


In ITj jil dz - - 1[* In 2 - £ , 


00 2 2m ~ 1 ( 2 2m - 1) B m 2m+l 


m^m+l ) ! 


<f>“ n,, ;i (A33) 


In (1-Z ) dZ = ^ (tt - <f> ) <j) + i [(f) (In 2 — 1 + In <f>) 


0 2m- 1 B o . 1 
2 m Zm+1 n 

2 ^To^ZTTT V J 


m 


i 1 m(2m+1)! 


(A34) 


Bernoulli's number is related to the Zeta function by the relation 


B 


n 


tt ( 2m ) ! 

0 2m-l 2m+l 

2 ir 


5 (2m) 


(A35) 


where c(2m) is the Zeta function given by 


s(2m) = k l 


1 

k 2m 


(A36) 


For m = 1, 2, 3, c(2m) = tt2/6, tt4/90, it6/945, respectively. With equations 
(A35) and (A36) 

m I I irsr+TJT «, 2m+1 = m -I , m( lm+7'J" (i) 2m+1 = k 4 / m I 1 SlW 

(lL.) 2n,+1 (A37) 

By partial fractions, l/m(2m+l) = 1/m -2/(2m+l), then the m-summation of 
equation (A37) is 



, also 


E 

k = 1 


1 n[l - (^) 2 ] = In sin <p - In 


thus the summation term of equation (A34) simplifies to 


2 2m ~ 1 B m ,2m+l 
2 .m(2m+l ) ! * 


m = 1 


lnsincf>d(j) + 4)“(j)ln4) 


In 4> - <p In sin $ + 2 tt“ (£■ - k tanh (A38) 


k = 1 


Similarly, the summation terms of equations (A31 ) , (A32), and (A33) are 


” (2 2m ' 1 - 1) B m 2m+l ,4 , 1n 2 

- f , m ( 2m+l ) ! ♦ * ♦ " * " k"= l" 


2tt e („ + k 


tan h-i ^ - 4k tanh-1 2 ^k) 


(A39) 


, the summation term of equation (A32) is the same as that of equation (A39) 
but with j - <t> replacing <f> 9 also 


E 

m = 1 


2 2m-l , 2 2m . , ) B m 2m+1 » r it . k n + li ) ln (i + 2±) + 

m(2m+l ) ! * = 2 k = 1 11 * k ¥k 


odd 


2cj), 


kd -tf) in 0 -&)] 


(A40) 


Example summations of trigonometric series. - With equations (A21), (A29), 
and (A30) 
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£ ^ n _2 = tanh _1 Z = [cos (b+2) <j> + i sin (b+2) <t>] 

n = 3 
odd 

( ^- tanh" 1 cos <p + i yjy ) 

then since the cosine series is the real part and the sine series the imaginery 
part 

y C0 ^ + b|<j) _ 1 cos (b+2) (j> tanh" 1 cos cf> - ^r£-r 
n = 3 n ”^ * 4 m 

odd 


~ = -- sin (b+2) cp tanh" 1 cos <p + 

odd 


where an alternative function for tanh” 1 cos <p is -In tan (<f>/2). 


With equations (A12), (A29), and (A30) 


cos (b+1) <j> + j sin (b+1) <p - j (cos b <p + i sin b 4 >) 


i sin 2 c)))(^- tanh" 1 cos <j> + i 4 j|y) 


,n+b 


n= E 3^ 

odd 


(1 - cos 2 


sin (b+2) cp (A41 ) 


cos (b+2) <p (A42) 


°° cos(n+b)(j) _ cos(b+l U sinfl r . 

„ 2 0 n(n-2) 2 2 L 

n = 3 

odd 


(b+1) <p tanh" 1 cos cp + 

(p 


cos (b+1) <j>] 


(A43) 
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I + Slgt [cos (b+1) * tanh-1 cos * - ^ 

odd 

sin (b+1) <f>] (A44) 


With equation (A28) divided by Z^ s and equations (A29) and (A30) 


,n-h 


n | 3 n(n-2)(n+lT 3 

odd 


z h z -k_ z k ! z Z ’* 5 + 

- 1 fi (Z l - ? - ?) In (1 - Z) g— 


(Z' 1 + \ - |) In (1 + Z) 


= j cos 2 - J sin 2 ( sln2 2 + ^ \ Sln ^On 2 sl ’n |- + i 
(|- - g | )] - j- cos j (cos 2 - i j sin <f> n 2 cos + 
i f] + 1 sin f 


then letting <j> = <j>/2, the trigonometric summations become 

« COs(2n-l)(f) COS d> ird> . q , d> . 0 , . o -io 

n l 3 n(n-2) (n+l ) = "3 ' ifiT sln * ' 3 Sln 3 * + sl " 2 * cos * ln 2 

odd 

sin <p - c ~ °y - ^ In 2 cos <j> (A45) 


* sin(2n-l )(f) sind> -n-cb ^ - o <j> ~ 9 

n | 3 n (n-2)(n+1) 3 ' 2|*f Sln * cos <f, - § cos 3 4 , + sin * cos 2 * 

odd 

In 2 cos <j> - — 3 ^ In 2 sin $ (A46) 
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Spanwise Loading Distribution 


This loading function is given in equation (21) and contains the summation 


2 f n sin n cf> (A47) 

n = 3 
odd 

Substituting equation (10) into equation (A47) results in 


E f n sin n 

n = 3 
odd 



odd 


sin(n-2)4>b sin nj> 
n(n-2) 


sin(n+2)4>b sin n 4> 


2 cos^b sin(n-l)4>b sin n <fi 2 cos^b sin(n+l)^b sin n <t> -» 
n(n-l) ~~ n(n+lT ■* 


(A48) 


The first term within the brackets, by trigonometric identities can be written 
as 

sin (n - 2 ) 4 ^ sin n <j> = ^ cos [n $ - (n - 2 ) 4 ^] - j cos [n <p + (n - 2 ) 4 ^] 

= y cos cos (n - 1 ) ip sin ip + sin (n - 1 ) ip -y cos i/j cos (n - 1 )ip + 
-j sin ip sin (n -1) (A49) 

where \p + and ip are defined as 

= 4)^ + 4), ip_ ~ 4> b - 4> (A50) 

Apply equation (A49) to equations (A43) and (A44) with b = -1, results in 
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n e 3 Sin ^n(n - ir 1n ^ = T [cos ' cos *- + I si " <* + - l*J> + 

odd 

sin 0 + sin ip_ (tanh” 1 cos ip_ - tanh” 1 cos ^ + )] (A51 ) 

The second term in equation (A48) is of the form obtained from equation 
(All). In trigonometric summation it becomes 

- cos(n+b)cj) _ cos(b-l)<t) cos(b+l)<t> sin<f> r n , + 

n | 3 n(n+2) 2 3 2 L 2]tf cos (b - 1 ) 4, + 

odd 

sin (b - 1) $ tanh" 1 cos ((>] (A52) 


co sin(n+b)cj> 

n 2 3 "<" +2 > ' 

odd 


sin(b-l )<j> 
2 


sin t b+D i . slnt t s . n (b . 


cos (b - 1) <f> tanh” 1 cos cj>] 


(A53) 


The sine product is expanded as before, but with b = 1, then the second term 
summation of equation (A48) becomes 


°° sin(n+2)<f>h sin n j> 

n = 3 
odd 


= \ [-cos + cos 4 j_ + j sin j (\fj + + $_) 


1? r 

sin j (^ + ~ */>_) - 2 sin " U [) + sin ip + sin ip_ ( tanh 1 cos + 
tanh” 1 cos ij> + )] (A54) 


The third term in equation (A48) is of the form of equation (A6), then 
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£ = cos (b + 1) <f> - ^yjy sin cos (b + 1) <f> + sin (b + 1) * 

n 3 
odd 

-cos (b + j) <p cos ^ In 2 cos + sin (b + ^-) <f> sin In 2 sin (A55) 

E S n(n-f)^ = sin ( b+ 1) <J> - 2 j|j sin | sin ( b + y) 4> - f cos (b + 1 ) <j> 
odd 

-sin (b + - 2 ~) 4> cos In 2 cos - cos (b + ^r) <j> sin In 2 sin ^ (A56) 

The sine product is expanded in terms of n - p then with b = in equations 
(A55) and (A56) the third term of equation (A48) becomes 


sintn-! . ) . ^ sin if M cos_*b [( ^ . . ++) sin ^ 


-2 cos<J> b 

■ 

E 

n = 3 

n(n-l) 


odd 

u - - 

, . *+ 

4 > + ) sin 2~ cos 

cos 2 



. *+ 
sin 

2 sin y~ 

sin ^ In 

4, 1 


cos 


cos 


2~ 


( A57 ) 


sin 


Similarly the fourth term in equation (A48) is developed as 

o a °° sin(n+l)fh sinn<i> _ cos ^b r o /, L , \ „ 1 / , 

2 cos4> b j A 2 ~ *-~ 2 Sln O sin 2 ^+ 


- * J - 


n = 3 
odd 


m|j_ 


" y_ y_ y+ r_ 

(y y -j- - - ip + ) sin y ~ cos 2 — + ('ll - ip_ - ij/ + ) sin cos ^— + 2 

ifr- ib+ 

ib. lb cos «— tJ j, \b sin 1F- 


cos 2“ cos 2“ 1 n 


cos 


2 

Tp- 

2~ K 

— 2 sin o— sin 75— In 

ifr+ 2 2 


sin 




(A58) 
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where $ and are defined in equation (A50). 

The desired summation of equation (A48) is obtained by adding together 
equations (A51), (A54), (A57), and (A58). Thus 

z f p sin n<f> = - s - n -y ^ ) (1+2 cos 2 <j> b ) si n^ + j (cos <j> b + cos §) 2 
odd 


tanh” 1 


sjnJb s in t + 1 , 

1 + cos cos 4> 2 ' 


4> b - cos 4> ) 2 


tanh 1 


sin 4>b sin <p 
1 - cos <j> b cos cp 


(A59) 


where from equations (2) and (6), n = cos <f» , n b = cos cf> b . Using hyperbolic 
function identities, equation (A59) in terms of n and n b becomes 


s f n sin n<f> = - i- (1 + 2 n h 2 ) (1 
n = 3 n J D 

odd 



1 +r U? n 

|n b +nj 



- n) z cosh" 1 ? - 

IVI 


(A60) 


The condition of wing bending moment taken about the midspan is n b = 0, then 
equation (A60) reduces to 


z f sin n<j> = - i (1 - n 2 )^ + n 2 cosh” 1 — (A61 ) 

n = 3 n J 

odd 

When either equation (A59) or (A60) is substituted into equation (21) the loading 
distribution becomes a closed function for any value of span station about 
which bending moment is taken. 
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The k-Factor Summation 


Combining equations (26) and (10) results in 


k 


oo 


n = 3 
odd 


" f n 


2 


00 f sin 2 (n-2)4>h sin 2 (n+2)4>h _ 2 sin (n-2)<l>bsin(n+2)4>b 

n | 3 1 n (n- 2 ) 2 n(n+ 2) 2 “ n(n- 2 )(n+ 2 ) 

odd 


4cosc{>k 


r sin(n-2)<f'hsin(n-l )$b sinfn+n<t>hsin(n+2)'t’h . si n(n-2)'t'hsin(n+l ) 
f- n(n-2T(n- T) n(n+2j(niTJ n(n-2l(n+TF 



2 sin(n -1 HhsinCn+l Hh + sin 2 (n -1 )*h + 
*"n(n-l ) (n+ 1 ) n(n-l ) 2 


sin 2 (n+l Hb -il (A62) 

n ( n+TP 

The summations are obtained by using equations (A22) through (A24) with the trig- 
nometric relations given in equations (A29) through (A34). An example the sixth 
term in equation (A62) has the form of the summation example of equation (A45). 

By trigonometric identity 


sin (n- 2 ) <j>. sin (n + 1 ) 


= 2 cos 


3 *b ■ 


cos ( 2 n - 1 ) 


(A63) 


then the desired summation is ^ cos 3$^ times 
minus times equation (A45) in which <p = <j>^. 


equation (A45) in which <p = 0, 
Thus 


°° sin(n- 2 ) 0 hsin(n+l )4>h 
E n( n- 2 ) (n+ 1 ) 


n = 3 
odd 


= sin 34 > k + j 2 S1 * n 


3 sin< 


cos 


sin 2 ^ cos^^ln sin^ + ^ cos 3 ^ In cos^ 


(A64 ) 


This mathematical procedure is used to develop the other summations of equation 
(A62). The results are listed as follows: 
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00 



odd 


sin 2 (n-2Hh 

n(n~2) 2 


= f ' \ sin 4(f, b^ ” \ sin2( l ) b ’ \ sin 2 tj) b cos 2 cj) b In 


tancj> b 


(A 65 ) 


n = 3 
odd 


sin 2 l 

! n+2 i 

tfh 

, nl 




; 't (<f, b ' 1 sin 4 <’b ) + l sin2,,, b 


4 >b “2" s i n 2 4>bCOS b In 


tan * b + Lr + - cos_.iib 


36 


18 


(A66) 


o “ sin(n-2)<^hsin(n+2)^h 
„ l 3 n (n-2)(n+2) " 

odd 




- 4 cos 2 <j) b ) sin 2 4> b + s i n 2 4> b cos b In 


tan<f>. 


(A 67 ) 




n = 3 
odd 

in2 


sin z 4> b cos 0 b In sin<}> b + cos 3 ci> b In coS(f) b 


(A68) 


2 

n = 3 
odd 


s i nl nj )tb?1n(n+2)* b = - ^ sin 3* b 


n(n-l ) (n+2) 

I 

In s i n 4) b ) + ^ cos 3 <f> k In cos<j> 


T 2 s i r ' 3< f’ b + 


\ sin 2 4> b coS(f) b ( 


(A 69 ) 


o “ sin(n-l )^hsin(n+l ) 4 >h 
’ n = 3 n(n-l ) (n+1 ) 

odd 


= sin 2 <j> b In sin<j> b - cos 2 (j) b In cos^ b (A 70 ) 
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co|rv) 



I 


~ sin*(n-1H h _ 

n = 3 "l"' 1 )* 
odd 


>b 1 1 

~2 4" ^5 ” g" n 2<j> b ” 2 ^ n2< ^b S1 " n( ^b ” 2^ 


cos 2 «f> b In cos<^ b 


(A71) 


00 sin 2 (n+l )^h 
n = 3 n^rHTj2 

odd 


‘ l +b + t sin 2<t, b ' 1 sin2<,, b 


-g sin 2 <t>b In sin<f> b - 2 


cos 2 $ b In cos^ 


(A72) 


Then by equation (A62) the k factor is the sum of equations (A65) through 
(A67), plus 4 cosc})^ times the sum of equations (A68) s (A64), and (A69), plus 

4 cos 2 ^^ times the sum of equations (A70) through (A72). Thus with an extensive 

reduction the k factor simplifies to 

2 

k = g (sin 6 cj> b - 3 sin 2 <j> b cos 2 ^ - 6 cos 4 ^ In cos 2 <t> b ) (A73 ) 

Since n b = coS(j) b> k can be written as 

k = -| (1 - 6 n b 2 + 3 n b 4 + 2 n b 6 - 6 n b 4 In n b 2 ) (A74 ) 


The functional behavior of k as n b -> 1 can be assessed by expanding the 
k factor in terms of 1 - n b - Then 


k = ^ (1 - V [1 - f O - n b ) + fo O * V 2 + 36 J 3 (n ~ 1 infe! h ) - ] 

odd 


(A75) 
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It can be seen that k decreases rapdily as n b approaches unity and remains 

positive. The summation terms are quite small. When an analytical summation 
of this series is made, equation (A75) converts back to equation (A74). 


Spanwise Center of Pressure, n c p 

The n of the spanwise loading distribution of equation (21) is given in 
cp 

equation (22) which contains the summation 


n-1 


-3 z 
n = 3 
odd 


Hi 


n 2 -4 


n-1 


y ~ 3 H) ^ 

L o n 2 -4 
n = 3 

odd 


sin(n-2)f h _ sin(n+2)tj>b _ 

n(n-2) ~ n(n+2) 


2 cosd>hsin(n-l )<f>h j. 2 cos<j)bsin.(n+]_) fbn 

RtHT] "WT) 

These summations are obtained by expanding the n 2 - 4 term in partial fractions, 
as 


-3 


1 H 

4 


- (- 
* [ r\-Z 



(A77) 


With equation (A77) the individual summations of equation (A76) appear as those 
of equations (A25) through (A27), which can be converted to trigonometric series 
by application of equations (A29) through (A34). The first summation by using 
equation (A25) with a = -2, b = c = 0, is 


n-1 


• z 

n = 3 
odd 


Hi 


sin(n-2)<f>h _ 
n(n-2) 2 


- Imag. (- \ tan 


n + 


4z 2 


tan l Z 


1 

2 


tan 1 Z 


dZ) = re- (2<j> k + sin 2* h ) - i sin* 


- 4 cos 


-<{>. tanh"* 1 s i n ^ (A78) 


The other summations are as follows: 
84 


J 



n-1 

« (-1 ) 2 sin(n-2)'f’h 

n f 3 n(n-2) (n+2) / 

odd 


j 2 s i ncf)^ - j sin^j^ cos 2 4> b + J sin^ cos : 



b cos 2^ b tanh" 1 


sin4> b 


(A79) 


n-1 

n ! 3 ^ L) iirJ]r^2f )il1 = - T? sin +b + I s1n *b cos2 i>b - f s1n *b cos3< f*b - 

odd 


^ cos 2 cf) b cos 24 > b tanh -1 s i n 


(A80) 


n-1 


-* ^ * § >«*> * i ><"»„ -> 2 - 

n = 3 
odd 


8 S1 ' n +b 


cos4> b - j cos 2 ^ tanh 1 sinc^ 


(A81) 


n 


n-1 

~ (-1 ) 2 sin(n-l )<f>h 

f 3 n(n-2)(n-1) 


odd 



(A82) 


n£L 

«? (-1 ) 2 sin (n-1 )<J>b 

n "f 3 n(n-l )(n+2) - 

odd 

j cos 3 4 > b tanh" 1 sin<f> b 


g- - j2 sin ^ b " f sin ^b C0s: 


+ ~ 


j sin^g cos<j> b + 


(A83) 
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n-1 


-£ 
n = 
odd 


(-1 ) 2 sin(n+] )4>h 
2 n(n-2)(n+l ) 



IT sin *b 


+ 


7T 

6 


s i n^f-, c °s Z( Pfo 



COSc}), 


+ y cos 2 ^ tanh" 1 sinfj 


(A84) 


n-1 

2 

f (-1 ) sin(n+l )4>h 
n = 3 n(n+l ) (n+2) 

odd 




COSc^ 


(A85) 


The f n summation term of equation (A76) is 3/4 times the sum of equations 
(A78) through (A81), plus 1.5 cos <j>^ times the sum of equations (A82) through 
(A85). Then, defining kj as follows: 


n-1 

2 

k i = - 3 l n 2- 4 f n = " 1 sin +b + if sin3 ‘f > b + 1 cos4<(, b tanh " 1 sin<f> b 
n - 3 

odd 

(A86) 

In terms of = cos^* k l is 

k i = i (1 • f n b 2) n b 2) ^ + 1 n b 4 cosh " 1 iT ( A87 > 


Induced Downwash Angle in the Wake 

The equation for induced angle is given in equation (23) which contains 
the summation 
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1 


sin * n = 3 
odd 


s nf sin nd> 
n 


... _J <*> r _ sin(n-2)4>b sin n<j» , sin(n+2)<j>h sin n± 

r i Mi L p_2 n-1-9 


sinc^ 


n = 3 
odd 


n+2 


2cos4 ) h sinfn-l)^ sin n<|> _ 2cos4>b sinfn+lKb sin n<H 
n-1 ~ n+1 J 


(A88) 


where f is defined in equation (10). By trigonometric identities the general 
sine product is formulated as 

sin (n + c)^ sin cos c^ cos mf> - ^ sin c^ sin n^_ - cos c^ ^ 

cos n^ + + j sin c<f> b sin m|< + (A89) 

where ip and are defined in equation ( A50 ) With c = -2 in equation (A89), 
the first summation of equation (A88) can be evaluated by using equations (A41 ) 

and (A42) in which b = 0, and 4> equals or i/> + . Then 

-z s1n l n z 2 J jb S1 ' . n n< P = _ II (i + . ) sin 2(f) - ~ (tanh _1 cos^ - tanh _1 

n = 3 n_Z 0 l^_l 

odd 

cosi|) + ) cos 2<f) (A90) 

Equation (A20) can be made into trigonmetric series by using the relations 
of equations (A29) and (A30). Thus 

v n+b i 

? -—^2 = - j cos (b + 1 )(f> - cos (b - 1)4> - sin (b - 2)<j> + j 

odd 

cos (b - 2)<f> tanh" 1 cos(f> +i [- y sin (b + 1 )<j> - sin (b - 1)4 + 

cos (b - 2)4> + sin (b - 2)<f> tanh” 1 coS(f>] (A91 ) 
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where the summations of cos (n + b)4> and sin (n + b)4> are given by the real 
part and the imaginary part, respectively, of equation (A91). With c = 2 in 
equation (A89) and using equation (A91), the second summation of equation 
(A88) becomes 


« sinCn+Z^b sin n<f> 

z g n+2 

n = 3 

odd 


3 Sln 


! 4> b sine), 


- - n 
8 u 


Vk 


-) sin 2 <f>+ 


(tanh -1 cos ip_ - tanh" 1 cosif/ + ) cos 2<p (A 92) 

In similar procedures, working with equations (16A), (17A), and (89A), the 
third and fourth summations of equation (88A) are evaluated as 


“ sin (n-1 )i > b sin n<f _ 
= o n-1 

odd 



sin<j> + | (1 


+ sine)) + ^ cosej) In 


sinto. 

T 

sinifi 


(A93) 


-Z s1n(n+l)»b s l n jlfc = s1n ^ cos^ sin* 

n = 3 ' ‘ ' 

odd 


sin<)) - 4- cos<p In 


s i mp + 
si mp 


1 

2 


4> b sine)) + ^ (1 



(A94) 


The required summation is the sum of equations (A90) and (A92), plus 
2cos<)> b times the sum of equations (A93) and (A94), then all divided by sincf>, 
thus 
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J 

sintf) 


n = 
odd 



sin n<j> 


2 

3 


(2 + cos 2 4> b ) sincj) b 


U b cos<j> b 




(cOSef) - COScf) b ) 



+ 


T^ )( * 



(A95) 


where fi is defined in equation (9). Shown in equation (A95) is a linear 
variation with spanwise coordinate n for n ^ n b » and is a constant with 
ti for n i n b - That is, the induced downwash in the wake is constant inboard 
of n b » then varies linearly with n at span stations outboard of n b - 


APPENDIX B 


EVALUATION OF I AND I * INTEGRALS 
n nn* 

Integral of Equation (67) 


The I integral is 




i = - 
n 7 t 


(cos^i -ti) cos n<t>i dfi 
(cos^-n) 2 + 4 2 


(Bl) 


This integral was originally evaluated, for arbitrary integer n, in reference 
14 in terms of a recurrence formula of a related integral. Here a simpler 
I recurrence formula will be developed. From trigonometric identities 

cos(n + 2)cj) 1 = 2 cos 2<f> x cos - cos(n - 2)<|> 1 
then integrating as in equation (Bl) 


(B2) 


n+2 


n t 


o 


cos 2^1 (cos^T-n) cos n<h d<f>i _ T 
(cos^i-n) 2 + £ 2 " i n-2 


(B3) 


In the integrant, cos 2<j>! can be divided by the denominator, and noting that 
2 cos4>icos n<|>i = cos (n + l)^ + cos (n - l)cj>i, then 
r 77 

9 

I 


^(cos^j-ri) cos n<f>i + 2 n(cos<f>i -n) [cos (n+1 )4 >t+ cos(n-l ) i J 
o l (cos<j>i-n) 2 + C 2 


n+2 ir 

( 1 +2 n 2+ 2 c 2 ) (cos^t-tOcos n<h 

{cos^y- n)2 + V 


Then the recurrence formula for I is 

n 


d *i - T n -2 


(B4) 


I „ = [2, n = 1] Rn, n = 0] + 4 n I . - 2(l+2 n 2 + 2 ? 2 ) I + 4„ 

n+d Lo, n t lJ Lo, n f oJ T n \ n (B5) 

n-l ‘ n-2 

Examination of equation (Bl) shows 


!_ n = r n » = l n^ for odd n ’ = "^(n) for even n (B6) 
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that is, I n is symmetric with n for odd n and antisymmetric for even n. Then 
from equation (B5) 


1 2 = - 2 n + 4 n I-, - (1 + 2 n 2 + 2 C 2 ) I Q 

1 3 = 2 + 4n - (3 + 4 n 2 + 4c 2 ) I-j + 4n Iq 
where based on the work in reference 14, I Q and 1^ are 

I = O-n) P~^-(l + n)r~^ 

0 [(pV )2 - A] h 


(B7) 

(B8) 


(B9) 


t _ -i (-n+n 2 + c 2 ) P ^“(n + n 2 + c 2 ) r - ^ 
1 1 ~ n i; r 

1 [(p 2 +r 2 ) 2 - 4]" 2 

C 2 (p~^+r~^) 

Xr 2 1 . 

C(p 2 +r^) - 4]^ 

where 


1 + 



(BIO) 


] 


(BID 


P = (1 - n) 2 + C 2 

r = (1 + n) 2 + c 2 = P + 4n 

Equations (B5) through (Bll) gives the evaluation of the I integral of 
equation (Bl) for all values of n, c, and n. Another method is, to define 

fTT 


J =1 
n 7 t 


cos n<h d4>i 

(cos*i-n)2 + ? 2 


(B12) 


From reference 14 


p + r 2 
E(p^-r^) 2 - 4] 3 * 


[(p^r 2 ) 2 - 


(Bl 3) 


(B14) 



Using a trigonometric identity, equation (Bl) can be written 


*rr 

[cos(n+1 ) ( j>i+cos(n-l )^i -2ncos njjJ dft] 
(cos^-nP- + c 2 

then 




2 (°n+l 


- 2ri J 


+ J ) 
n-1 


(Bl 5) 


Since 


J_ n = J n i J n (-n) = - J n U) for odd n; J r/' n ^ = J n ^ for even n ( B16 ) 

that is, J n is antisymmetric with n for odd n and symmetric for even n. 

Then 

W J 0 (B17) 

Equation (B2) can be used to establish a recurrence formula for J n in a manner 
similar to that for I . Then 


4, n = 01 + 4n J , -2(lf2r, 2 +2 c 2 )Jf4 n J , -J ? 
0, n f OJ n 1 n 

then 



J 2 = 2 + 4 n 0, - (1 + 2 n 2 + 2c 2 ) J Q 

= 4n ~ (3 + 4n 2 + 4c 2 ) J-j + 4n Jq 
Combining equation (B15) for n = 1, and equation (B19) gives 


(Bl 8) 


(Bl 9) 
(B20) 


I 1 = 1 + nJ 1 - (n 2 + e) J Q = 1 + nl 0 - C 2 J Q (B21) 

Inserting equations (B13) and (B14) into equations ( B 1 7 ) and (B21) leads to 
the I Q and I-j values of equations (B9) and (BIO). The J^ recurrence formula 
of equation ( Bl 8 ) , together with equations (Bll), (B13), and ( Bl 4 ) provide 
an evaluation of the J n type of integral (eq. B12) for arbitrary n» C s and n. 
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With J n known, a second method, by applying (B15), is available for evaluat- 
ing I n . 


Integral of Equation (70) 


The I nn * integral is 


I 


nn* 


i r 

— I n * s i n<f> sin n«j> d<f> 


(B22) 


Jj o 

where I * is that in equation (Bl) with n = n*. Examination of equation 
(B22), taking into account the symmetric and antisymmetric characteristics 
of I * in equation (B6), shows 


*0n* O’ *-n,n* *nn*’ *n,-n* *nn*’ 

*n odd, n* even ^n even, n* odd 

which means that I nn * has a value only when n and n* 
n* are both even integers. 

A recurrence formula for I nn * can be made by making the integration of 
equation (B22) on equation (B5). Then 


= 0 


(B23 ) 


are both odd, or n and 


; 1 

n,n*+2 it 


n 

~2 , n* = 1' 

1 [4 cos cj>. 

n* = 0] 

1 

0 

_0, n* f 1. 

J " Lo 

n* f Oj 


n*+l 


COS d) - 


2(1 + 2^ 2 + 2 cos 2 <|)) I n * + 4 1^* ^ coscj> - 1^* sin«|> sin n<j> d<f> 


(B24) 


Inserting the trigonometric identities 

2 sin n4> cos<j) = sin (n + 1)4> + sin (n - l)<f> 

4 sin n<j> cos 2 <j) = 2 sin n$ + sin (n + 2)<j> + sin (n - 2)^ 

into equation (B24), then the integrations can be made in terms of I-values with 

lower n* integers, which is a recurrence formula. The recurrence formula for 

I * is 
nn* 
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'l , n = 1 , n* = 1 ' 
n ,n*+2 [o, n / 1 , n* / 1_ 


‘1, n = 2, n* = 0] ?T 

_0, n f 2, n* 1 oj *n+1,n*+l 


21 


i i,i ~ 4(1 +5^) I . - I 0 + - I 0 * + 21 * , + 21 

n-l,n*+l v ’ nn* n+2,n* n-2,n* n+1 ,n*-l n-l,n*-l 


" I n,n*-2 (B25) 

The I nn * terms for low values of n and n* are obtained from equation (B25) by 
taking into account the identity relations given in equation (B23). Thus 
letting n* be zero in equation (B25) results in 


l n2 " 2 _0, n,t 2_ + 2I n+l ,l +2I n-l ,r 2 ^ 1+?2 ^ I n0'2 ^+2,0 2 1 n-2,0 


*02 ! 1 2 J n odd, 2 ~ 0 


! 22 ' 2 + 2I 3 1 +2I ir 2 ( 1+ ^) I 20T ! 40 

*42 = 2I 51 +2I 31" 2 ^ 1+c 2 ^ I 40"2 ^OT *20 
Letting n* be one in equation (B25) results in 


> (B26) 


I 

I 

I 

I 


n3 

03 

13 

33 


*1 , n = 1 “ 
.0, n f 1 


2 I n+l,2 +2I n-l,2' (5+4?2)I nr I n+2,T I n-2,l +2I n+l,0 +2 ^n-1,0 


^00 ^23 *n even, 3 ” ^ 

1 + 2 1 22 " 4 (l+c 2 )I-|-j - I 31 + 2I 20 

2I 42 + 2I 22 " 3] " 1 5 i " l -[-[ + 2I 40 + 2I 2Q 


>(B27) 


J 

For higher values of n*, equation (B25) is used directly, without the bracketed 
terms since n* is greater than one. 
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With the recurrence formula of equation (B25) the I nn * terms can be 
evaluated to very large values of n*. However, as can be seen in equations 
(B26) and (B27), the recurrence formula starts with a low value of n*, 
namely n* = 0 and 1. Needed are evaluations of I q and 1^ . 

For n « 1* by binomial theory, the radicals of equations (B13) and (B14) 
can be expanded into a power series of n» then 


J o = 


|c|(l+e 2 y 


TT [1 + (A - C 2 ) ( 


_n — \2 + t± 

1 +C 2 ' '8 


3c 2 + C 4 ) (i^) 4 + 


% - ir ?2 + ™ ?z - ?6) ( ^ )6] 


l+c 


(B28) 


- 1 r 1 ■ 


J, = 1 r- libr + (j - 2; 2 ) 


) 3 + (#- 


§ C 2 + 3 ? 4 ) (^) 5 + 


(fe - ^ ?2+ 15 s 4 -4c 6 ) (i^) 7 ] 


15 

2 


(B29) 


The influence coefficients I Q and I-j are obtained by combining equations (B28) 
and (B29) with equations (B17) and (B21). Thus 


= 


13 i 

2 C r n 

r = 1 r 
odd 


(B30) 


where 


C-, = 


Co 


- c 


( 1^ 2 ) 3/2 

-Lcl 


(in 2 ) 


7/2 v 2 


(I- 


- c 2 ) , 


c c = 


c 7 = 


Hi 


(l+c 2 ) 

.. zlsl 


(l+c 2 ) 


1 

(11 . 

5c 2 

11/2 

I 

,35 
16 " 

105 

15/2 

8 ' 


(B31) 


2,21 4 

C 2 + o - C 4 


c 6 ) 

J 


and Cg> C-ji* and 3 are determined from the simultaneous solution of three 
equations, at n = . 8 , .9, and .95 of 


7 

n 9 Cg + n 1 1 c-j 1 + n 1 3 c-j 3 = Ig(n) - £ C r n ] 

odd 


( B32) 
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where I Q (n) is given in equation (B9). The inverse of equation (B32) for 
n = .8, .9, and .95, is 


29241 A V’ 8 ^ 

23104 A V* 9 ^ 

82944 aI 0 ( - 95 > 

1785 

,8 9 

629 

.9 9 

3885 

T95^ 

13700 

Alg(.8) 

, 61700 

Al 0 (.9) 

46400 

AI q ( .95) 

357 

.8 y 

629 

“.99 

777 

. 95 s 

CO 

o 
1 — 1 

< 

o 

o 

o 

CO 

II 

40000 A V' 9 ^ 

, 32000 Al 0 ( ' 95 > 

357 

.8 9 

629 

.9 y 

' 777 

.9P~ 


(B33) 


where 


AI 0 (n) = I q ( n ) - 2 c r n r 

r = 1 
odd 


(B34) 


And for I 


12 i 

T 1 = * C r n 

r = 0 r 
even 


(B35) 


where 


5 0 = 1 “ 


(H C 2 ) 

-3kl 
2 2(l+5 2 ) 5/2 


1/2 


Co = 


- C 


H = 


(He 2 ) 


(He 2 ) 


,15 5 

9/2 ^8 "2 c ^ 


/ 3 5 35 9,7 (. \ 

(it - — e 2 + 2 c 4 ) 


TTJY V 16 “ 4 


(B36) 


29241 A M' 8 ^ 23104 A V' 9 ^ 82944 A M’ 95 ^ 

1785 ~7P 


629 


79^ 


3885 795^ 


_ 1 3700 Al i ( * 8 ) _ 61700 Al i ( - 9 ) 46400 Al l^’ 95 ^ 

c 10 357 “P 


629 


. 9 8 


777 T9P" 


_ _ 8000 a M* 8 ) 40000 Al l^' 9 ^ , 32000 A M' 95 ^ 

c 12 357 .8 8 " 629 .9 8 777 .95 8 


(B37) 
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where 


AMn) = Un) - e c 
r = 0 1 

even 


(B38) 


For evaluating I nQ and I pl , equations (B30) and (B35) are inserted into 
equation (B22). Since n = cos<p, and sincf> sin n<j> = 0.5 [cos(n-l)<j> - cos 
(n+l)<|>], then 


t - J_ 

nO 2 77 E . c r 
r = 1 

odd 

which integrates to 


cos r 4> [cos(n-l)<f> - cos(n+l)$] d<f, 


I = ^JL 

nO 4 


[..rj .f* 


n= 


2 ’ 2 

1 4 

4’ ^ 




n=* 


I6> 2 

1 4 
4 ’ q 

16’ 6 




7 n= 
32* 2 


32 


, 4 


1_ c 

32’ D 

M* 8 


JLa 

4 


21 

128 s 

n= 

2 

+ ^-L 
+ 4 

33 

256 s 

— 

n= 

2 

3 



165 


16 s 

4 


1024 s 

4 

27 

a 


55 

a 

256 s 

O 


512 s 

0 

1 

32 s 

8 


11 

256 s 

8 

1 

256 s 

10 


5 

512 s 

10 




1 

12 


1 

1024 s 




429 

4096’ 

143 

1024 s 

429 

4096 s 

13 

256 s 
65 

4096 s 
3 

1024 s 
1 

4096 s 


n= 

2 

4 

6 

8 

10 

12 

14 


(B39) 


where the brackets indicate that a value exists at the cited n, but is zero 
for any other n integer. The c r 's for odd r are given in equations (B31) and 
(B33) . 

Similarly, with equations (B35) and (B22) 
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I = lo_ 
nl 4 


[«.n 




£jb. 

4 


7_ 

64* 

7_ 

32’ 

5_ 

32 s 

7 

128’ 

1 

128 s 


n= 

1 

3 

5 

7 

9 


1 r 

2 s 1 

1 3 

2 s o 


+ S_L£L 
4 


2] 

256 s 
45 
256’ 
75 

512 s 

35 

51 2 * 
9 

51 2 * 
1 

51 2 ’ 


n= 

1 

3 

5 

7 

9 

11 


4 - 

4 


r 


1 

4 s 


n~ 

1 


+ ll2_ 
4 


— 3 
8 s 

1 5 
8 s 


r 33_ 

51 2 * 
297 
2048 ’ 
275 
2048 * 
77 

1024* 

27 

1024 s 
11 

2048 s 

1 

2048 s 




n= 

1 

3 

5 

7 

9 

11 

13 


5 n=: 

12’ 1 

12 s 3 
5 

32 s 5 

jL 7 

32 s ' 


(B40) 


where c r 's for even r are given in equations (B36) and (B37). 

Mlg ncal exam 2 le_gf I nn * for - An application of the equations 
for evaluating I nn * leads to the following values: From equations (B31) and (B36) 


-.357771 


C 3 = -.286217 
C 5 = -.100748 


Cq = .552786 
C 2 = -.429325 
C 4 = -.228973 


C 7 = .042497 

t 6 = -.025645 

From equations (B9) and (B10) 


I 0 (.8) = -.442322 

q ( .8) = .199389 

I Q ( - 9) = -.529048 

I 1 ( . 9) = .100090 

I Q ( .95) = -.569635 

I,(.95) = .051389 

Using the c's of equation (B41) 



(B41) 


(B42) 
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n 

7 r 

E C n 

r - 1 r 

odd 

n 

S r 

r = 0 r 
even 

.8 

-.456861 

.8 

.177508 

.9 

-.569810 

.9 

.041175 

.95 

-.633557 

.95 

-.040031 


(B43) 


By equations (B34) and (B38) the aIq(ti) and Al-j(n) terms are given by the 
differences between the values of equations (B42) and (B43). Then with 
equations (B33) and (B37) 


ii 

cr> 

kj' 

.075149 

e 8 = .051299 




11 

.107226 

C 1Q = .191347 


> 

(B44) 

13 = 

-.086553 

C ]2 = -.105817 

J 




With the s values established in equations (B41) and (B44) the 
and I^i coefficients are determined from equations (B39) and (B40). These 
coefficients are listed in the first two columns of table 12. In table 12 
the I nn * for n* > 1 are determined from the recurrence formula of equations 
(B25) through (B27). 
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APPENDIX C 


EVALUATION OF P AND P * INTEGRALS 
n nn* 


Integral of Equation (125) 

The recurrence formula of equation (125) is obtained in a similar manner 
as done for I in appendix B. With equation (B2), equation (125) becomes 


P 


n +2 


2 cos 2 ^ (cosycos 4>i-n) cos njj d 4 > t _ p 

7 T ( cos i -n cosy ) 2 + n 2 sin^y ~ n -2 

o 


(Cl) 


In the integrant, cos 2<j>! can be divided by the denominator, and noting that 
2 cos <}> 2 cos rhp 1 = cos(n+l )</>!+ cos ( n-1 ) 4 >i , equation (Cl) becomes 


P 


n+2 


-4n, n = 0 
2 cosy, n = 1 
0 , n > 1 


n-1 


which is the recurrence formula for P . Examination of equation (125) shows 


P_ n = P n » P n (-n) = p n (n) for odd n ’ p n (~n) = “ P n^ for even n ( C3 ) 


Then from equation (C2) 

?2 ~ + 4^ cosy P-j - (1 + 2r\ 2 ) Pq (C4) 

Pg = 2 COSy + 4f) COSy Pg “ + 4r | 2 ) P-j + 4iy COSy Pg (C5) 

where 

Pg = cosy - n Jg (C6) 

P-j = COSy + n COS2y J-, - n 2 COSy Jg (C7 ) 


where Jg and J^ are given in equations (B13) and (B14), respectively, but with 
the parameters ncosy Cor n and nsiny for 5 in p and r of equation (Bll). Then 

p = 1 - 2 ncosy + n 2 1 

> (C 8 ) 

r = 1 + 2 ncosy + n 2 
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Equations (C 6 ) and (C 7) are derived from a parameter similarity with the Jq 
and J-j integrals. Thus 


p _ cosy 

0 “ 7T 


fir 


cos<h djj 

(cos^i-ncosy ) 2 + n 2 sin 2 y 


f IT 


d^i 

(cosc^-ncosy ) 2 + n 2 sin 2 y 


which comparing with the J n integral of equation (B12) results in 


Pq = cosy J-| (ncosy, nsiny) - n Jq (ncosy, nsiny) 


that is, these J’s are functions of ncosy and nsiny rather than the n and r, 
of equation (B12). Equation (C7) is derived in a similar fashion. 

In terms of the ncosy, nsiny parameters the J n integral is 


,7T 


J n 


cos n4> i d4 > 7 

(cos^i-ncosy ) 2 + n 2 sin 2 y 


(CIO) 


Using the equation (B2) relation, a recurrence formula can be developed, thus 


J 


n +2 


[4, n - 0" 
LO, n f Oj 


+ 4n cosy J n+1 - 2(1 + 2n 2 ) J p + 4n cosy J n _-| - J n _ 2 

(cn ) 


Since J_ n = J n , then 


J 2 = 2 + 4 n cosy J ] - (1 + 2 n 2 ) J Q (C12) 

Jg = 4n cosy J 2 - (3 + 4n 2 ) J-| + 4 n cosy Jq ( Cl 3 ) 

where J Q and J-j are given in equations (B13) and (B14) in which p and r are 
those in equation (C 8 ). 


When n and 5 are replaced by ncosy and nsiny in equations (B18) through 
(B20) , then equations (Cl 1 ) through (C13) are reproduced or checked. Equation 
(125) can be written as 


P 

n 


fTT 


^cosx [ CO s(n+l )$]_ + cos(n-l) 4 > 1 ] - ncos n<f>]^ d^j 
(costjij-ncosy ) 2 + n 2 sin 2 y 
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P p = \ ( J n+1 c °sy - 2n J n + J n _ 1 cosy) (C14) 

Equation (C14) provides a second method (to eq. C2) for evaluating P . 

Integral of Equation (128) 

The recurrence formula for P nn * is made by making the integration of 
equation (128) on equation (C2) in which n is replaced by n*. Then in a 
derivation similar to that leading to equation (B25), the recurrence formula 
is 


P * , o ~ fcosy, n = 1 , n* = ll +2 cosy (P + P ) - 

n, n*+2 i n = p n* = n n +l , n*+l n-1, n*+l 

L -1 ’ n £ ’ n u J 4P * - P x . * 

nn n+2, n* 

-P « * + 2 cosy ( P , i *i + P i * -, ) - P 

n-2, n* n+l, n*-l n-1, n*-l n, n*-2 

(Cl 5) 

where the brackets indicate that this term is zero for values of n and n* 
other than those listed. The relations of equation (B23) apply also to P pn *- 
Then from equation (Cl 5) 

P n2 = b n=2 ^ + 2 COSy ^ P n+1 ,1 +P n-1 ,1 ^ " 2P nO ' 2 P n+2,0 ' 1 P n-2,0 

P 02 = P 12 = P n=odd,2 = 0 * ^ C16 ^ 

P 22 = " 2 + 2 COSy ^ P 31 + P 1 1 ^ " 2P 20 " 2 P 40 


P n3 “ £ C0ST ’ n 1] + 2 cos T( p n+1 j2 +P n-l ,2 ) *“ 5P nr P n+2,l" P n-2,l +2 C0Sy (P n+l ,0 +P n-l ,0^ 


Pm = Po-3 = p , Q = 0 

03 23 n=even,3 


P ^2 = cosy + 2 cosy P ^ “ ^ P H “ P 31 + 2 C0S Y P 20 


(Cl 7 ) 


Values of P n Q and are needed to start this recurrence formula. From 
equations (128), (C6), and (C7) 
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nr 


_ 1 


no 


P 00 P n=odd,0 


(cosy J i - J q cos<f>) sin<j) sin n < p dtf> 
= 0 


(C18) 


where the zero values result because the term inside the parenthesis is anti- 
symmetric about 4> = tt/2, so when sin n<j> is symmetric, the integral is zero. 


P 

P 


nl 

01 


= [|T cosy, n=l ] 

= P n = 0 

n=even,0 


+ 


fir 


1_ 

2ir 

J 0 


(cos 2y J^-cosy JqCOS<j>) sin 2cf) sin n<f> d<f> 

(Cl 9) 


where as before the term inside the parenthesis is antisymmetric. These inte- 
grals can be evaluated by the expansion method used in appendix B, or by a 
quadrature formula given in reference 15, which has 

f 17 

= jjjjp [-J- f{* 0 ) + \ f(<t> M+1 ) + K fU )] (C20) 

. 0 w" 1 

where M is an odd integer and $ = yir/(M+l). The integrant of equation (Cl 8) 

is symmetric about (j> = tt/ 2 when n is even, also it is zero at <f> = tt/ 2. Then 
with equation (C20), equation (C18) becomes 


n=even ^ 
P n0 = M+T 


M-l 

2 

2 (cosy J-j - Jq cosef)^ ) sin<£ sin n<f> 


(C21 ) 


where Jq and J-j of equations (B13) and (B14) are evaluated with the p and r of 
equation (C8) at n = cos4> . 

The integrant in equation ( Cl 9 ) is symmetric when n is odd, also it is 
zero at <f> = tt/2. Then with equation (C20), equation (C19) becomes 


n=odd i 1 

P nl = [j cosy, n=l] + jjfpf 


M-l 

2 

z (cos 2y J, - cosy J n cos<t> ) sin 2<f> sin n<j> 
y=] 1 u y y y 

(C22) 


where as before, ^ = utt/(M+ 1), and Jq and J-j are evaluated with n = cos<f>^ 

in equation (C8). 
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TABLE 1. - CONSTANTS OF PLANAR WING 
MINIMIZATION SOLUTION 


% 

k 

ki 

fi 

1 1 — 1 

0 

2/9 

1/6 

4/3 

1/6 

.05 

.218893 

.165424 

1.181253 

.185306 

.10 

.209570 

.161760 

1.039157 

.201673 

.15 

.195126 

.155839 

.907010 

.215131 

.20 

.176851 

.147886 

.784747 

.225361 

.25 

.156042 

.138175 

.672280 

.232109 

.30 

.133952 

.127012 

.569493 

.235213 

.40 

.090328 

.101 679 

.392358 

.230218 

.5 

.053025 

.074704 

.251841 

.210548 

.6 

.025900 

.048928 

.145912 

.177501 

.7 

.009611 

.026977 

.071919 

.133638 

.8 

.002196 

.010978 

.026398 

.083197 

.9 

.000157 

.002159 

.004719 

.033267 

1.0 

0 

0 

0 

0 




TABLE 2. 


SPANWISE FUNCTION, k Q , OF PLANAR WING 
MINIMIZATION SOLUTION 



0 

.2 

• 4 

.6 

.8 

.9 

.96 

1 

0 

-1/3 

-.239745 

-.059810 

.161042 

.406024 

.534748 

.613552 

2/3 

.05 

-.325360 

-.238882 

-.060532 

.159478 

.403914 

.532429 

.611122 

.664168 

.10 

-.308364 

-.236111 

-.062646 

.154815 

.397606 

.525492 

.603852 

.656692 

.15 

-.286239 

-.230808 

-.065990 

.147144 

.387163 

.513998 

.591800 

.644294 

.2 

-.261029 

-.221316 

-.070279 

.136621 

.372713 

.498050 

.575063 

.627069 

.3 

-.206572 

-.183739 

-.079716 

.107996 

.332454 

.453424 

.528138 

.578723 

.4 

-.152579 

-.139061 

-.083345 

.071782 

.278779 

.393352 

.464721 

.513248 

.5 

-.103773 

-.096063 

-.066443 

. 032844 

.214672 

.320444 

.387268 

.433013 

.6 

-.063166 

-.059116 

-.044010 

.001076 

.144553 

.238432 

.299239 

.341333 

.8 

-.012386 

-.011761 

-.009496 

-.003586 

.020040 

.071386 

.113089 

.144000 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 
















o 

CO 


TABLE 3. - FC 
OF 


\ b 

n 

0 

.05 

.10 

3 

2.6667 

2.6500 

2.6005 

5 

-.3810 

-.3710 

-.3418 

7 

.1270 

.1199 

.0996 

9 

-.0577 

-.0515 

-.0366 

n 

.0311 

.0266 

.0147 

13 

-.0186 

-.0149 

-.0053 

15 

.0121 

.0089 

.0011 

17 

-.0083 

-.0055 


19 

.0059 

.0035 

-.0018 

21 

-.0044 

-.0022 

.0016 

23 

.0033 

.0014 

-.0022 

25 

-.0026 

-.0008 


27 

.0020 

.0005 

-.0018 

29 

-.0016 

-.0002 

.0016 

31 

.0013 

0 

-.0013 

33 

-.0011 

.0001 

.0011 

35 

.0009 

-.0002 

-.0009 

37 

-.0008 

.0002 

.0007 

39 

.0007 

-.0002 

-.0005 

41 

-.0006 

.0003 

.0003 

43 

.0005 

-.0003 

-.0002 

45 

-.0004 

.0003 

.0001 

47 

.0004 

-.0003 

0 

49 

-.0003 

.0003 

-.0001 

51 

.0003 

-.0002 

.0001 


CJ O') oo 
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TABLE 4. - AERODYNAMIC CHARACTERISTICS OF WING WITH 
WINGTIP ZERO-SLOPE LOADING 


% 

t = t ! 



ttA 

2C l 

a w 

"1 

C b xlO 

— 

e 

C b 2 

— 

V /e 

C bc 


n cp 


e X|U 


S 

eq. (29) 

o<n<n b 


n=l 

0 

3/2 

.31831 

3 



-1.71239 

.79577 

2/3 

.94989 

.84375 

4/3 

.05 

1.50564 

.31871 

2.77854 



-1 .71506 

.67774 

.66835 

.68727 

.77779 

1.38697 

.10 

1.52278 

.31987 

2.58241 



-1.72315 

.57298 

.67296 

.43785 

.71342 

1.44322 

.15 

1.55209 

.32176 

2.40776 



-1 .73687 

.48077 

.68025 

.33979 

.65224 

1.50128 

.20 

1.59472 

.32432 

2.25145 



-1.75652 

.40005 

.68977 

.23202 

.59496 

1 . 561 01 

.25 

1.65247 

.32751 

2.11092 



-1 .78262 

.32979 

.70121 

.15510 

.54193 

1.62220 

.30 

1.72794 

.33127 

1.98405 



-1.81589 

.26900 

.71431 

.10130 

.49323 

1.68473 

.40 

1.94837 

.34033 

1 .76446 



-1 .90813 

.17218 

.74466 

.03981 

.40838 

1.81340 

.5 

2.30940 

.35119 

1.58160 



-2.04600 

.10296 

.77954 

.01360 

.33859 

1.94643 

.6 

2.92969 

.36358 

1 .42748 



-2.25408 

.05573 

.81813 

.00380 

.28158 

2.08343 

.7 

4.11847 

.37726 

1.29620 



-2.58537 

.02573 

.85983 

.00077 

.23511 

2.22412 

.8 

6.94444 

.39206 

1.18332 



-3.18000 

.00887 

.90424 

.00009 

.19722 

2.36830 

.9 

18.11177 

.40782 

1 . 08547 



-4.60451 

.00149 

.95102 

.00000 

.16608 

2.51588 

1.0 

00 

. 

.42441 

1 



— 00 

0 

1 

0 

.14063 

2.66667 







TABLE 5. - AERODYNAMIC CHARACTERISTICS OF WING WITH 
(C b c /C b> =1-1 LOADING 


n b 

t 

eq. (46) 

\p 

ttA 

2C l a w 

C b xl0 

e 

C h 2 

& y i n 2 

C b 2/ e 

e 

< c b 2/e >c 

o*n*n b 

n=l 

0 

6/11 

.38583 

1.72727 

.01368 

.96458 

.93798 

.99192 

.88109 

.1 

.45077 

.39347 

1.46843 

.19389 

.75176 

.9591 6 

.58921 

.86164 

.2 

.40339 

.39909 

1 .31656 

.30272 

.56771 

.97203 

.33157 

.85023 

.3 

.38650 

.40358 

1.22011 

.37016 

.41199 

.98038 

.17313 

.84298 

.4 

.39488 

.40737 

1.15494 

.41060 

.28384 

.98611 

.08170 

.83809 

.6 

.51216 

.41378 

1.07473 

.43113 

.10556 

.99325 

.01122 

.83206 

.8 

1.09270 

.41932 

1.02885 

.34228 

. 01 91 0 

.99738 

. 00037 

.82861 

1.0 

CO 

.42441 

1 

— CO 

0 

0 

0 

.82645 






TABLE 6. - AERODYNAMIC CHARACTERISTICS OF WING WITH 
(C bc / C b) = 1.2 LOADING 


1 

t 

eq. (46) 

n cp 

_7L! 

2C 

A 

— a 
L W 

C b xlO 

e 

C 2 

b xl fl2 

C b Ve 

CKn^ 

i n=i 

1 

« X 1 u 

e 

0 

1 

.35368 

2.33333 

-.80826 

.88419 

.81818 

.95553 

.84877 

.] 

.82642 

.36768 

1 .85878 

-.47787 

.68911 

.87479 

.54284 

.79384 


.73955 

.37800 

1.58036 

-.27834 

.52040 

.91180 

.29701 

.76162 


.70858 

.38622 

1 .40353 

-.15471 

.37766 

.93698 

.15222 

.74115 


.72395 

.39317 

1.28405 

-.08568 

.26019 

.95480 

.07090 

.72732 


.93896 

.40491 

1.13701 

-.04293 

.09676 

.97768 

.00958 I 

.71030 

.8 

2.00328 

.41508 

1.05288 

-.20581 

.01751 

.99126 

.00031 

.70057 

1.0 

00 

.42441 

1 

— 00 

0 

1 

0 

.69444 














TABLE 7. - AERODYNAMIC CHARACTERISTICS OF WING WITH 
(C bc / C b ) = V3 LOADING 


n b 

t 

eq. (46) 

n cp 

irA 

TTTT- a 

2C l w 

C b xlO 

e 

C h 2 

b xi n 2 

C b 2 /e 

e xl ° 

^b 575 ^ 

o<n<n b 

n=l 

0 

3/2 

.31831 

3 

-1.71239 

.79577 

2/3 

.94989 

.84375 

.05 

1.34912 

.32969 

2.59365 

-1 .43281 

.70501 

.71510 

.69506 

.78661 

.10 

1.23963 

.33931 

2.28817 

-1.21680 

.62020 

.75641 

.50852 

.74365 

.15 

1.16208 

.34755 

2.05402 

-1 .04915 

.54133 

.79145 

.37026 

.71072 

.20 

1.10933 

.35479 

1.87054 

- .91751 

.46836 

.82126 

.26710 

.68492 

.25 

1.07708 

.36125 

1.72410 

- .81371 

.40124 

.84672 

.19014 

.66433 

.30 

1.06287 

.36712 

1.60530 

- .73207 

.33989 

.86857 

.13301 

.64762 

.40 

1 .08593 

.37755 

1.42607 

- .62085 

.23417 

.90374 

. 06068 

.62242 

.5 

1.18738 

.38677 

1.29903 

- .56610 

.15031 

.93044 

.02428 

.60455 

.6 

1.40844 

.39517 

1.20551 

- .56439 

.08709 

.95113 

.00797 

.59140 

.7 

1.87073 

.40299 

1.13454 

- .62858 

.04292 

.96746 

. 001 90 

.58142 

.8 

3.00492 

.41041 

1.07932 

- .80872 

.01 575 

.98055 

.00025 

.57366 

.9 

7.51495 

.41753 

1.03546 

-1.32543 

.00282 

.99121 

.00001 

.56749 

1 .0 

<30 

.42441 

1 

— oo 

0 

1 

0 

.56250 




TABLE 8. - SPANWISE LOADING DISTRIBUTION, c„c/C,c . OF WING 

X, L dv 

WITH SPECIFIED CONDITION OF LOADING* 





Wingtip zero-slope loading 




n b 

n 

t 

0 

.2 

.4 

.6 

.8 

.9 

.96 

0 

3/2 

1 . 90986 

1.69614 

1 .27164 

.77254 

.29867 

.10982 

.02840 

.1 

1.52278 

1 .87112 

1.69605 

1 .27827 

.77846 

.30140 

.11088 

.02869 

.2 

1.59472 

1 .80325 

1.68781 

1.29773 

.79667 

.30988 

.11419 

.02957 

.4 

1 . 94837 

1.65175 

1.58552 

1.35644 

.87613 

.34900 

.12965 

.03371 

.6 

2.92969 

1.50886 

1 .46357 

1 .31740 

1.01538 

.44042 

.16731 

.04397 

.8 

6.94444 

1.38276 

1.34940 

1 .24390 

1.04396 

.65763 

.27986 

.07653 

elliptic load- 








ing 

either t=0 

1.27324 

1 .24751 

1 .16694 

1 .01859 

.76394 

.55499 

.35651 

or 

j 












(C bc /C b ) = 4/3 loading 




0 

3/2 

1.90986 

1 .69614 

1 .27164 

.77254 

.29867 

.10982 

.02840 

.1 

1.23963 

1.75994 

1 .61265 

1.25757 

.82311 

.38741 

.19346 

.08964 

.2 

1 .10933 

1 .64193 

1.55380 

1.25792 

.86422 

.44808 

.24836 

.12908 

.4 

1 . 08593 

1 .48420 

1.43590 

1.27256 

.9391 9 

.53267 

.31793 

.17659 

.6 

1 .40844 

1 .38651 

1.35138 

1.23928 

1 .01705 

.60841 

.36862 

.20625 

.8 

3.00492 

1 .32063 

1 .29160 

1.20024 

1 .02957 

.71794 

.43594 

.23536 





( C bc /C b) 

.2 loading 




0 

1 

1.69765 

1.54660 

1 .23674 

.85456 

.45376 

.25821 

.13777 

.1 

.82642 

1.59771 

1 .49094 

1 .22736 

.88827 

.51292 

.31397 

.17860 

.2 

.73955 

1 .51903 

1 .45170 

1.22760 

.91568 

.55337 

.35057 

.20489 

.3 

.70858 

1 .45961 

1.40993 

1 .23286 

.94065 

.58398 

.37668 

.22309 

.4 

.72395 

1 .41388 

1 .37311 

1.23735 

.96566 

.60976 

.39695 

.23657 

.6 

.93896 

1.34876 

1 .31676 

1 : 21 51 6 

1 .01756 

.66025 

.43074 

.25634 

.8 

2.00328 

1 .30483 

1 .27691 

1 .18914 

1 .02591 

.73327 

.47562 

.27574 





< C bc' C b> ■ 

1.1 loading 




0 

6/11 

1.50474 

1 .41065 

1.20501 

.92912 

.59475 

.39311 

.23720 

.1 

.45077 

1 .45022 

1.38029 

1.19990 

.94751 

.62702 

.42353 

.25947 

.2 

.40339 

1.40731 

1.35889 

1.20003 

.96246 

.64908 

.44349 

.27381 

.3 

.38650 

1.37489 

1.33611 

1.20290 

.97608 

.66578 

.45773 

.28373 

.4 

.39488 

1.34995 

1.31602 

1.20535 

.98972 

.67985 

.46879 

.291 08 

.6 

.51216 

1.31443 

1.28529 

1.19325 

1 .01803 

.70739 

.48722 

.30187 

.8 

1.09270 

1.29047 

1.26335 

1.17905 

1.02258 

.74722 

.51170 

.31245 


* Equation (63) and table 9 apply when comparing with elliptic loading with same 
values. 


113 



TABLE 9. - COMPARISON OF AERODYNAMIC CHARACTERISTICS DUE 
TO OPTIMIZED SPANWISE LOADING WITH THAT DUE TO 
ELLIPTIC LOADING OF EQUAL LIFT AND BENDING MOMENT 
FOR BENDING MOMENTS AT THE SAME n h = n hr NORMAL- 
IZED SPAN STATION D 


n b 

n bc 

b/be 

C bc /C b 

t 

eq. (59) 

n cp 

a /a 

w' wc 

C b xlO 

e 

D i 


n=l 

D ic 

0 

.95 

-.31579 

.44675 

.64149 

1 .74075 

1 . 1 1 688 

.97832 

1.13259 


1.00 

0 

.42441 

1 

1 

1 .06103 

1 

1 


1.05 

.28571 

.40420 

1.25256 

.43842 

1 .01051 

.98218 

.92348 


1.10 

.54545 

.38583 

1 .42750 

.01131 

.96458 

.93798 

.88109 


1.15 

.78261 

.36905 

1 .54516 

-.31392 

.92264 

.88020 

.85906 


1.20 

1.00000 

.35368 

1.62037 

-.56129 

.88419 

.81818 

.84877 


1 .25 

1 .20000 

.33953 

1 .66400 

-.74874 


.75758 

.84480 

* 

4/3 

3/2 

.31831 

1 .68750 

-.96322 

.79577 

2/3 

.84375 


1.60 

2.25000 

.26526 

1.56250 

-1.19867 

.66315 

.47059 

.83008 


2.00 

3 

.21221 

1.25000 

-1 .10619 


.33333 

.75000 

.1 

.95 

-.26097 

.43656 

.80754 

1 .62515 


.98593 

1.12385 


1 .00 

0 

.42441 

1 

1 

.82693 

1 

1 


1.05 

.23612 

.40820 

1.12958 

.52404 

.78756 

. 98845 

.91763 


1.10 

.45077 

.39347 

1 .21358 

.16024 

.75176 

.95916 

.86164 


1.15 

.64676 

.38001 

1 .26434 

-.11841 

.71907 

.91 940 

.82243 


1.20 

.82642 

.36768 

1 .29082 

-.33185 

.68911 

.87479 

.79384 

* 

1.44322 

1.52278 

.31987 

1.23982 

-.82729 

.57298 

.67296 

.71342 


1.60 

1 .85945 

.29676 

1.14541 

-.90828 

. 51 683 

.57985 

.67367 


2.00 

2.47926 

.25420 

.89409 

-.85840 

.41347 

.43703 

.57204 

.2 

.95 

-.23354 

.43907 

.90496 

1.55533 

.65735 

.99045 

1.11872 


1.00 

0 

.42441 

1 

1 

. 62448 

1 

1 


1.05 

.21130 

.41115 

1.05743 

.57575 

.59474 

.99217 

.91419 


1.10 

.40339 

.39909 

1.08807 

.25018 

.56771 

.97203 

.85023 


1.15 

.57878 

.38809 

1.09958 

-.00033 


. 94407 

.80094 


1.20 

.73955 

.37800 

1.09747 

-.19329 


.91180 

.76162 


1.35 

1.15042 

.35221 

1 .04405 

-.54240 


.81034 

.67712 

* 

1.56101 

1.59472 

.32432 

.92396 

-.72084 


.68977 

.59496 


2.00 

2.21866 

.28516 

.68527 

-.70875 

.31224 

.53460 

.46764 

.4 

.95 

-.22862 

.43428 

1 .00864 

1.48613 

.32866 

.99530 

1 .11326 


1.00 

0 

.42441 

1 

1 

.31223 

1 

1 


1.05 

.20684 

.41549 

.98064 

.62700 

.29736 

.99615 

.91053 


1.10 

.39488 

.40737 

.95450 

.33934 

. 28384 

.98611 

.83809 


1.15 

.56657 

.39996 

.92423 

.11670 

.27150 

.97182 

.77807 


1.20 

.72395 

.39317 

.89170 

-.05950 

fKmEl 

.95480 

.72732 


1 .40 

1 .24106 

.37086 

.75864 

-.43490 

wEBm 

.87787 

.58119 

* 

1.81340 

1 .94837 

.34033 

.53657 

-.58026 

.17218 

.74466 

.40838 


2.00 

2.17185 

.33069 

.46304 

-.56043 

.15611 

.70123 

.35652 

.6 

.95 

-.29651 

.43057 

1.06009 

1 .47296 


.99773 

1 .11056 


1.00 

0 

.42441 

1 

1 

.11611 

1 

1 


1.05 

.26827 

.41884 

.94253 

.63675 

.11058 

.99814 

.90872 


1.10 

.51216 

.41378 

.88821 

.35631 

.10556 

.99325 

.83206 


1.15 

.73484 

.40915 

.83722 

.13898 

mmm 

.98621 

.76672 


1 .20 

.93896 

.40491 

.78959 

-.02981 


.97768 

.71030 


1.60 

2.11266 

.38054 

. 51 1 04 

-.52601 

.07257 

.89638 

.43578 
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TABLE 9. - Continued. 


n b 

b/b c 

t 

n cp 

a /a 

w' wc 

C b xl0 

e 

D. 

l 

^bc 

c b</ c b 

eq. (59) 

Oin^ri^ 

n=l 


D 1c 

* 

2.08343 

2.92969 

.36358 

.32886 

-.51929 

.05573 

.81813 

.28158 

.8 

.95 

-.63261 

.42736 

1.08953 

1.52995 

.02211 

.99912 

1 .10706 


1.00 

0 

.42441 

1 

1 

.02101 

1 

1 


1.05 

.57236 

.42175 

.92073 

.59454 

.02001 

.99928 

.90768 


1.10 

1.09270 

.41932 

.85029 

.28288 

.01910 

.99738 

.82861 


1.15 

1.56778 

.41711 

.78744 

.04259 

.01827 

.99463 

.76023 


1.20 

2.00328 

.41508 

.73117 

-.14292 

.01751 

.99126 

.70057 


1.60 

4.50737 

.40341 

.48358 

-.66917 

.01313 

.95729 

.40806 

* 

2.36830 

6.94444 

.39206 

.21097 

-.56696 

. 00887 

.90424 

.1 9722 


* Values for wingtip zero-slope loading, t of equation (29). 


115 



TABLE 10. - COMPARISON OF AERODYNAMIC CHARACTERISTICS DUE TO OPTIMIZED SPANWISE 
LOADING WITH THAT DUE TO ELLIPTIC LOADING OF EQUAL LIFT AND BENDING 
MOMENT FOR BENDING MOMENTS AT THE SAME y b = y b SPAN STATION 


n b 

u 

JD 
\ II 
jQ 

t 

"cP 

“•/“wc 

C b xlO 


D i 







D ic 

b cV /b 

C bc /C b 

eq. (60) 


oin^n b 

n=l 


0 


Same values as 

those in 

table 9 for n b = 0. 



.1 

.95 

-.33009 

.44708 

.72796 

1 .76210 

.88198 

.97767 

1 .13334 


1 .00 

0 

.42441 

1 

1 

.82693 

1 

1 


1.05 

.29820 

.40394 

1.18810 

.42334 

.77720 

.98170 

.92393 


1.10 

.56887 

.38536 

1.31499 

-.01429 

.73206 

.93649 

.88250 


1.15 

.81558 

.36842 

1 .39699 

-.34668 

.69092 

.87765 

.86155 


1.20 

1 .04134 

.35292 

1.44592 

-.59876 

.65327 

.81483 

.85226 

* 

1.32306 

1.52278 

.31987 

1.47525 

-.98439 

.57298 

.67296 

.84889 


1 .6 

2.32882 

.26453 

1.33594 

-1.23617 

.43856 

.46804 

.83460 


2.0 

3.08624 

.21253 

1.05177 

-1 .12976 

.31224 

.33376 

.74903 

.2 

.95 

-.37441 

.44791 

.78247 

1 .82513 

.67717 

.97581 

1.13550 


1.00 

0 

.42441 

1 

1 

.62448 

I 

1 


1.05 

.33663 

.40328 

1.14664 

.37926 

.57711 

.98035 

.92521 


1.10 

.64065 

.38420 

1 .24194 

-.08874 

.53432 

.93233 

.88643 


1.15 

.91632 

.36690 

1 .29987 

-.44149 

.49552 

.87071 

.86842 


1.20 

1.16718 

.35116 

1.33051 

-.70659 

.46022 

.80585 

.86175 

* 

1.29721 

1.59472 

.32432 

1 .33796 

-1.04383 

.40005 

.68977 

.86154 


1.6 

2.56054 

.26370 

1.17554 

-1 .33827 

.26413 

.46307 

.84356 


2.0 

3.32803 

.21553 

.90292 

-1.18815 

.15611 

.33798 

.73969 

.4 

.95 

-.60668 

.45059 

.84428 

2.11139 

.35584 

.96782 

1.14487 


1.00 

0 

.42441 

1 

1 

.31223 

1 

1 


1 .05 

.53345 

.40139 

1.09687 

.18483 

.27388 

.97494 

.93034 


1.10 

1.00378 

.38110 

1 .15193 

-.41177 

.24008 

.91658 

.90166 


1.15 

1 .41938 

.36316 

1.17725 

-.84580 

.21020 

.84604 

.89375 


1 .20 

1.78725 

.34729 

1.18142 

-1.15809 

.18376 

.77608 

.89481 


1.22396 

1.94837 

.34033 

1.17782 

-1.27372 

.17218 

.74466 

.89641 


1.6 

3.56422 

.27060 

.93689 

-1 .68748 

.05603 

.46566 

.83887 


2.0 

4.19758 

.24327 

.66174 

-1 .31633 

.01050 

.38587 

.64789 


2.5 

4.34371 

.23697 

.43269 

-.87735 

0 

.36978 

.43269 

.6 

.95 

-1.44623 

.45445 

.87421 

2.88793 

.14592 

.94861 

1 .16806 


1.00 

0 

.42441 

1 

1 

.11611 

1 

1 


1.05 

1 .20278 

.39944 

1.06621 

-.30472 

.09133 

.96388 

.94101 


1.10 

2.19802 

.37877 

1.09150 

-1 .19123 

.07081 

.88879 

.92986 

* 

1.14423 

2.92969 

.36358 

1.09029 

-1 .72164 

.05573 

.81813 

.93358 


1.15 

3.01594 

.36179 

1 .08889 

-1.77684 

.05396 

.80933 

.93428 


1.20 

3.68188 

.34796 

1 .06752 

-2.14552 

.04023 

.74013 

.93827 


1.40 

5.21515 

.31612 

.89845 

-2.44520 

.00863 

.58671 

.86960 


5/3 

5.63377 

.30742 

.65593 

-1.89273 

0 

.54883 

.65594 

.8 

.95 

-7.85059 

.46099 

.87840 

6.34397 

.03472 

.88079 

1.25800 


1 .00 

0 

.42441 

1 

1 

.02101 

1 

1 


1.05 

5.43663 

.39908 

1.03720 

-2.06115 

.01150 

.93905 

.96590 

★ 

1.06860 

6.94444 

.39206 

1.03627 

-2.78482 

. 00887 

.90424 

.96847 


1.10 

8.94496 

.38274 

1.02159 

-3.62327 

.00537 

.85055 

.97166 


1.15 

10.94770 

.37341 

.97467 

-4.22658 

.00187 

.79164 

.95516 


1.25 

12.01966 

.36841 

.84307 

-3.99032 

0 

.75915 

.84305 


* Values for wingtip zero-slope loading, t of equation (29). 


116 




CM "3- CO 00 


TABLE 11. - SPANWISE LOADING DISTRIBUTION*, c c/C. c , 

a L av 

FOR CONDITIONS OF TABLE 10 


% 

Nv n 
t 

0 

.2 

.4 

.6 

.8 

.9 

.96 




(b/b c ) = 1.20 loading 




.1 

1.04134 

1.68209 

1.55424 

1 .24307 

.85438 

.44764 

.25129 

.13233 

.2 

1.16718 

1 .66115 

1.56977 

1 .26267 

.85617 

.43161 

.23237 

.11722 

.4 

1.78725 

1.62045 

1.55757 

1 .34077 

.88791 

.38331 

.16482 

.06040 

.6 

3.68188 

1.56936 

1.51905 

1.35603 

1 .01456 

.35735 

.06778 

-.03628 




(b/b c ) = 1.15 loading 




.1 

.81558 

1.59345 

1 .48775 

1.22657 

.88998 

.51621 

.31713 

.18093 

.2 

.91632 

1.57778 

1 .50051 

1.24209 

.891 08 

.50304 

.30171 

.16865 

.4 

1 .41938 

1.54898 

1.49375 

1.30499 

.91481 

.46166 

.24513 

.12135 

.6 

3.01594 

1 .51580 

1.46993 

1.32183 

1 .01529 

.43089 

.15590 

.03476 

.8 

10.94770 

1 .44589 

1 .40814 

1.28826 

1.05858 

.59634 

.12126 

- . 08487 




(b/b c ) = 1.10 loading 




.1 

.56887 

1 .49659 

1 .41508 

1.20853 

. 92888 

.59115 

.38908 

.23404 

.2 

.64065 

1.48616 

1 .42440 

1 .21948 

.92944 

.58153 

.37791 

.22516 

.4 

1.00378 

1 .46824 

1 .42165 

1.26457 

.94520 

.55017 

.33586 

.19020 

.6 

2.19802 

1 .45002 

1 .40961 

1.27983 

1 .01618 

.52122 

.26413 

.12202 

.8 

8.94496 

1 .41430 

1.37876 

1.26607 

1.05126 

.62700 

.20060 

-.00413 


(b/b c ) =1.05 loading 


.29820 

1.39032 

1 .33535 

1.18874 

.97157 

.67337 

.46802 

.29231 

.33663 

1 .38512 

1.34046 

1 .19455 

.97175 

.66809 

.46194 

.28749 

.53345 

1.37687 

1 .34006 

1.21883 

.97959 

.65033 

.43854 

.26813 

1.20278 

1.36997 

1 .33622 

1.22871 

1 .01727 

.63112 

.39583 

.22819 

5.43663 

1 .35898 

1 .32728 

1 .22719 

1 .03845 

.68071 

.33960 



.13732 


* The loadings 
of table 10, 


for n b = 0, and 
are the same as 


for wingtip zero-slope loading with b/b c 
those in table 8. 
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TABLE 12. - I nn * s INDUCED DRAG INFLUENCE COEFFICIENT FOR s=l/2, FROM APPENDIX B 


0 

1 

2 

3 

0 

0 

0 

0 

0 

.211037 

0 

-.023182 

-.126495 

0 

.107717 

0 

0 

-.069557 

0 

.058029 

-.017041 

0 

-.039669 

0 

0 

-.002698 

0 

-.022936 

.002017 

0 

.001243 

0 

0 

.001782 

0 

.001958 

.000806 

0 

.001031 

0 

0 

.000244 

0 

.000516 

-.000008 

0 

.000026 

0 

0 

-.000049 

0 

-.000025 

-.000037 

0 

-.000023 

0 

0 

-.000013 

0 

.000006 

-.000005 

0 

.000006 

0 

0 

0 

0 

.000020 

0 

0 

.000003 

0 

0 

0 

0 

.000005 

6 

7 

8 

9 

0 

0 

0 

0 

0 

-.000042 

0 

-.000382 

.000418 

0 

-.000639 

0 

0 

.000550 

0 

-.000640 

-.008969 

0 

.000549 

0 

0 

-.005641 

0 

.000423 

.010127 

0 

-.003581 

0 

0 

.005505 

0 

-.002419 

-.004694 

0 

.002362 

0 

0 

-.003203 

0 

.000485 

.000900 

0 

-.002706 

0 

0 

.000371 

0 

-.002310 

-.000118 

0 

-.000040 

0 

0 

-.000337 

0 

.000077 

-.000190 

0 

-.000506 

0 

0 

-.000391 

0 

- 


0 

0 

-.019769 

0 


.032101 

0 

-.013384 

0 


.001766 

0 

.000243 

0 


-.000037 

0 

.000022 

0 


.000031 

0 


5 


0 

-.000442 

0 

-.013710 

0 

.017975 

0 

-.007837 

0 

.001378 

0 

.000066 

0 

-.000079 

0 

.000013 

0 

.000036 






















































TABLE 13 . - n *- D_* INDUCED DRAG INFLUENCE COEFFICIENTS OF PLANAR-WING 
it nn 

WINGLET WITH $ Q = 5 tt / 32 , FOR EQUATIONS ( 172 ) AND ( 175 ), (DATA FROM TABLE 

II OF REF. 9) 


Y 

N^n* 

1 

3 

5 

7 

9 

11 

deg 

n 







90 

1 

.03335 

.00666 

-.00167 

-.01243 

-.06192 

-.10534 


3 

-.00333 

-.17809 

-.16927 

-.05646 

-.02895 

-.06658 


5 

-.00258 

-.16334 

-.09138 

.08085 

.11343 

.03307 


7 

.03630 

.02441 

.11900 

.18434 

.11512 

-.00789 


9 

.03156 

.12523 

.17174 

.08391 

-.03334 

-.07625 


11 

-.02442 

.06277 

.04171 

-.11282 

-.16338 

-.02769 

75 

1 

.02254 

.01464 

.00887 

-.00752 

-.04600 

-.06860 


3 

-.02833 

-.11801 

-.09743 

-.03368 

-.02510 

-.04331 


5 

-.02362 

-.10743 

-.05579 

.05572 

.06810 

.01018 


7 

.02184 

.01 044 

.07264 

.11569 

.06939 

-.01221 


9 

.01822 

.08039 

.10960 

.04817 

-.02089 

-.03748 


11 

-.01661 

.04548 

.02751 

-.07520 

-.09798 

-.00506 
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TABLE 14. - T n , WING ROOT BENDING MOMENT/LIFT 
INFLUENCE COEFFICIENTS OF EQUATION (171) 


n 







cb 

0 

1 

3 

5 

7 

9 

11 

5tt/32 

.075139 

.184765 

.198455 

.120903 

.009690 

-.066188 

Tr/4 

.285398 

1/2 

1/6 

-1/6 

-1/10 

1/10 
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